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1. Introduction 

The physics program at the Large Hadron Colhder (LHC) rehes heavily on the ever 
increasing theoretical control over modeling high-energy proton collisions. The detailed 
theoretical understanding not only increases the reach in new physics and particle 
searches, but also allows to study the fundamental dynamics and properties of particles. 
Formulating new observables for addressing specific physics questions is a typical task 
which relies on quantitative reliable theoretical input. In the long run, high statistics 
measurements combined with precision prediction from theory will allow systematic 
probes of fundamental particle theory at ever deeper levels. 

The theory of Quantum Chromodynamics (QCD) concisely describes the collisions 
of protons, however, the dominant dynamics differ depending on observables and the 
regions of phase-space [1, 2, 3, 4, 5]. Here we have in mind proton collisions with 
large momentum transfer that are typical for the production of heavy particles. These 
include the Higgs, top, new particles within theories of supersymmetry as well as 
more conventional Standard Model processes at large scattering angles. To a good 
approximation such scattering processes factorize into the long-distance dynamics of 
quarks and gluons within protons, short distance hard interactions between these 
partons and, finally, the formation of hadrons and observable jets from the emerging 
partonic states and remaining proton fragments. Monte Carlo event generators deal with 
all aspects of the multi-layered simulation of the proton collisions. For some purposes, 
i.e. sufficiently inclusive observables, it is accurate to use a simpler, purely partonic 
description of events. To this end one combines final-state partons into observable jets, 
consistently ignoring corrections from showering and hadronization. Numerical methods 
are commonly used for the evaluation of differential cross sections, being well suited for 
comparisons to experimental data. 

The hard scattering process is a central stage in the simulation of proton collisions. 
It is described through scattering amplitudes which are accessible through first principle 
computations in quantum field theory. The complexity of scattering process allows only 
a perturbative approach, based on the expansion in the strong coupling asifJ^ji). ^ ^^^^ 
step towards a theoretical understanding of QCD is the evaluation of cross sections 
at leading order (LO) in the strong coupling as{l.i%)- Many tools [6, 7, 8, 9, 10] 
are available to generate predictions at leading order. Some of the methods applied 
incorporate higher-multiplicity leading-order matrix elements into parton-showering 
programs [11, 12], using matching (or merging) procedures [13, 14]. 

The truncation of the perturbative expansion introduces an explicit dependence 
on the unphysical renormalization scale leading to a theoretical uncertainty. QCD 
cross sections can have strong sensitivity on higher-order corrections, motivating the 
challenging quest for perturbative corrections. Next-to-leading (NLO) order predictions 
significantly reduce renormalization- and factorization-scale dependence - a feature that 
becomes more significant with increasing jet multiplicity (see e.g. [15]). In addition, NLO 
corrections take into account further physics effects including initial state radiation. 
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parton merging into jets and a more complete account of partonic production channels. 
Fixed-order results at NLO can also be matched to parton showers [16] with the prospect 
of complete event generation at next-to-leading-order in the strong coupling. The value 
of a first principle understanding of scattering processes in addition to the increased 
quantitative control motivates the quest for cross sections at NLO and sometimes beyond 
this [17, 18]. 

In this contribution we will describe loop-level on-shell [19, 20, 21] and unitarity 
methods [22, 23]. Our main focus will be on generalized unitarity [24, 25, 26] and its 
extensions for the numerical computation [27, 28, 21, 29] of hard scattering matrix 
elements. In addition, we will discuss numerical on-shell loop-level recursions [21]. 
The central strategy of these approaches is to make maximal use of universal physical 
principles and mathematical structures in order to describe complex multi-particle 
processes including quantum corrections. Numerical algorithms based on unitarity and 
on-shell methods are efficient, numerically stable and can be automated for generic 
scattering processes. The rapid recent developments in this field have already led 
to many new studies of complex scattering processes at the LHC [30, 31, 32, 33, 
34, 35, 36, 37, 38]. Several further implementations of unitarity methods have been 
reported [39, 40, 41, 38, 42]. For related analytic approaches see the chapter in this 
review by Britto [43] . 

For processes at the LHC with complicated final states, computations can be very 
challenging. In principle, one can use the traditional Feynman diagram representation 
of amplitudes for numerical evaluation. However, even at leading order, the number of 
Feynman diagrams grows more than factorial as the number of final-states increases. 
Computation times scale accordingly and refined approaches are needed. 

At leading order recursive approaches allow to reduce the growth in complexity. 
Two central strategies are commonly used for numerical computations: Off-shell 
recursion relations [44, 45, 46], based on Dyson- Schwinger equations, optimise the 
reuse of recurring groups of Feynman graphs. In contrast, on-shell recursions [47, 48] 
take advantage of the remarkable simplicity of the physical scattering amplitudes 
(see e.g. [49]). The simplicity arises in part from symmetry properties of tree 
amplitudes [50, 51, 52] that are present in QCD-like theories; see the chapter in this 
review by Brandhuber, Spence and Travaglini [53]. For most practical purposes the 
efficiency of the two approaches is comparable, depending on the explicit realization of 
the algorithms. 

At next-to-leading-order additional challenges arise, in particular, for the virtual 
corrections due to the loop-momentum integration. NLO cross sections are built from 
several ingredients: virtual corrections, computed from the interference of tree-level 
and one-loop amplitudes; real-emission corrections; and a mechanism for isolating and 
integrating the infrared singularities in the latter. Automated approaches [54, 55, 56] to 
deal with these issues are available based on subtraction methods [57, 58, 59]. Recursive 
methods are effective for computations of real-emission corrections. Such methods, 
however, are not directly applicable to virtual corrections. Traditional methods evaluate 
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the loop integrals of Feynman diagrams (sec e.g. [60, 61, 62]), and have to overcome 
two central challenges: growth of the number of Feynman diagram expressions and the 
evaluation of tensorial loop integrals, while maintaining gauge invariance. Means to 
deal with tensor integral reductions [63, 64] as well as strategies to recycle substructures 
have been shown to reduce complexity inherent in Feynman diagram approaches [65, 61]. 
For a more complete discussion of important NLO computations along these lines see 
ref. [66]. 

The unitarity method [22], in contrast, constructs loop amplitudes from on-shell 
tree amplitudes; gauge invariance is built in and maintained throughout computations. 
In addition, complexity arising from large numbers of Feynman diagrams is avoided by 
recursive methods for tree evaluations. Similarly, loop-level recursions [19, 20, 67, 21] 
construct amplitudes efficiently using purely on-shell lower-point input. Effective 
numerically stable implementations of these on-shell methods have been demonstrated 
by various groups [21, 68, 35, 40, 39, 41, 38, 42]. Beyond this, unitarity approaches have 
already been applied to state-of-the-art NLO computations [36, 30, 31, 32, 69, 34, 37]. 

In more detail, numerical unitarity based approaches use a combination of 
methods. Scattering amplitudes are naturally split into two parts; a part with 
logarithmic dependence on kinematic invariants and a rational remainder. Typically, 
unitarity approaches in strictly four dimensions are used for the computation of the 
logarithmic parts, although, on-shell recursions [20] may as well be applied in certain 
cases. At present, there are three main choices for computing the rational remainder 
within a process-nonspecific numerical program: on-shell recursion [21], /^-dimensional 
unitarity [23, 29], and a refined Feynman-diagram approach [70, 38]. We will discuss 
here numerical unitarity approaches in four and D dimensions following refs. [21, 29] as 
well as numerical loop-level on-shell methods [21]. 

Several recent developments allow us to use a purely numerical approach at loop 
level. A key tool is generalized unitarity [24, 25] which imposes multiple unitarity 
cuts (more then two propagators) and gives a refined system of consistency relations 
that is easier to solve. In addition, generalized cuts allow for a hierarchical approach; 
computing coefficients of four-point integrals, in turn three-point and, finally, two-point 
etc. integrals from successively cutting four, three and two propagators. (For the related 
maximal cut technique for multi-loop amplitudes sec [71].) Analogous approaches may 
be applied in D dimensions [23] (see also [72]), where higher point integrals have 
to be considered. Further simplifications arise from working directly with the loop 
integrand as opposed to the integrated loop amplitude. The unitarity method then 
turns into a purely algebraic approach. The starting point for this approach is a generic 
representation of the loop integrand [27], whose free, process dependent parameters 
are to be determined by unitarity relations. A particularly useful parametrisation of 
the loop integrand was given by [27] (see also ref. [28]) and extended in [29] to D 
dimensions. Importantly, despite the restrictions imposed by on-shell conditions in 
unitarity cuts, sufficient freedom remains in loop-momentum parametrization in order 
to uniquely determine all integral coefficients for renormalizable theories and beyond. 
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Modern on-shell and unitarity methods may be set up to take advantage of 
refined physical properties and formal structures of scattering amplitudes. We will 
discuss the uses of structures like spinor helicity [73, 74], analyticity properties, color 
decomposition [75, 76, 77] and supersymmetry properties [78, 24] in order to make 
computations efficient. The scattering amplitudes are then decomposed into a fine set 
of gauge invariant pieces (primitive amplitudes) , which are computed individually and 
eventually assembled into the full matrix element. This approach leads to excellent 
numerical stability and can be further exploited [32] for caching and efficiency gains 
through importance sampling (as used in color-expansions). In addition to aiming for 
efficiency it can be helpful to use methods which are easy to automate within existing 
frameworks [70, 38] or fulfill further computational constraints [40]. 

Furthermore, a numerical approach to amplitudes requires attention to numerical 
instabilities induced by round-off error. A natural way to deal with round-off errors is 
to setup a rescue system which monitors numerical precision and invokes an alternative 
computational approach when checks fail. A convenient rescue strategy [79, 21] for 
unitarity based approaches is the use of higher precision arithmetic. The advantage of 
a fine split-up of loop amplitudes into gauge invariant subparts is that one can setup a 
very targeted and thus efficient rescue system. 

The present chapter of this review is organized as follows. In section 2 we 
discuss a representative example for next-to-leading-order multijet computations at 
hadron colliders pointing out the importance of on-shell and unitarity methods. In 
section 3 we discuss the key properties of in matrix elements that can be exploited for 
the computation of loop amplitudes. Finally, numerical unitarity approaches will be 
discussed in section 4 and loop- level recursions in section 5. We end with conclusions 
and an outlook. 



2. NLO Predictions for Hadron Colliders 

As an example, to point out the key features of NLO QCD predictions we will focus on 
processes of massive vector-boson production in association with jets. In particular, we 
focus on recent progress due to the use of numerical unitarity approaches. 

The production of massive vector bosons in association with jets at hadron colliders 
has been the subject of theoretical studies for over three decades [80, 81, 82, 83, 84, 85, 
86, 87]. These theoretical studies have for example played an important role in the 
discovery of the top quark [88]. The one-loop matrix elements for W + 2-jet and Z + 2- 
jet production were determined [24] via the unitarity method [22] (see also ref. [89]), 
and incorporated into the parton-level MCFM [90] program. Studies of W production 
in association with heavy quarks have also been performed [91, 92, 93, 94, 95, 96]. 

Beyond this, the numerical unitarity approach allowed to include additional final- 
state objects. Studies of iy-(-3-jet production can be found in [30, 31, 32, 69] and Z + 3- 
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jet in [34]. The state-of-the-art in perturbative QCD for hadron coUiders arc currently 
parton-lcvcl next-to-leading order computations with up to five final-state objects. The 
first and only such process to be computed so far is + 4-jet production [36]. More 
generally, several important QCD processes with four final-state objects have been 
computed to date [97, 35, 98, 37]. 

Processes of Z- and 1^-boson production in association with jets have a particularly 
rich phenomenology at the electro-weak symmetry breaking scale, being important 
backgrounds to many searches for new physics and particles, for Higgs physics, and 
will continue to be important to precision top-quark measurements. Decays of the 
massive vector-bosons into neutrinos mimic missing energy signals and are of particular 
importance for supersymmetry searches [99, 100, 101]. The clean leptonic decays of 
the vector-bosons open a high resolution view of underlying QCD dynamics. Inclusive 
production cross section provides valuable information about parton distribution 
functions as well as fundamental parameters of the Standard Model. The signal 
of vector-boson production in association with jets per se includes physics of jet- 
production ratios [102] including comparative studies of W, Z and photon production. 
Experimental studies of vector-boson -|- jet production at the Tevatron were published 
by the CDF and DO collaborations [103, 104, 105] as well as at the LHC by the ATLAS 
collaboration [106]. 

2.1. Validation & Prediction 

Before turning to the theoretical and technical issues it is useful to assess how good 
the results are by comparing to experimental data. The quantitative impact of NLO 
corrections can be directly vahdated against data from Tevatron pp coUisions. Fig. 1 
compares the Et distribution of the third-most energetic jet in CDF data [103] to the 
NLO and LO predictions for W + 3-jet production [30, 32]. (For a similar analysis 
using a leading color approximation see also [31].) The upper panels of fig. 1 show the 
distribution itself, while the lower panels show the ratio of the LO value and of the 
data to the NLO result. The NLO predictions match the data very well, and uniformly 
with well matching slope. The central values of the LO predictions, in contrast, have 
different shapes from the data. The change of shape between LO and NLO is attributed 
to an imprecise scale in the coupling constant, that governs the production of the softest 
observed jet, in the leading order computation which gets corrected once loop effects 
are included as discussed in refs. [107, 31] . 

Scale-dependence bands indicate rough estimates of the theoretical error. Those 
are obtained by varying the renormalization- and factorization scales by factors of two 
around a central scale. For such scale variations, the dependence is of the order of ±40% 
for W + 3-jet processes at LO. At NLO, the scale dependence shrinks to ±10%, and 
we obtain a quantitatively reliable answer. A more detailed discussion can be found in 
refs. [30, 32]. 

As another example, we consider predictions for the LHC. For the inclusive 



Susy Theories and QCD: Numerical Approaches. 



7 




Figure 1. The measured Et distribution of the softest observed jet in inclusive 
W + 3-jet production, compared to the NLO prediction [30, 32]. In the upper panels 
the NLO distribution is the solid (black) histogram, and CDF data points are the 
(red) points, whose inner and outer error bars, respectively, denote the statistical and 
total uncertainties (excluding the luminosity error) on the measurements. The LO 
predictions are shown as dashed (blue) lines. The lower panel shows the distribution 
normalized to the full NLO prediction. The scale-dependence bands are shaded (gray) 
for NLO and cross-hatched (brown) for LO. 

Reprinted fig. 1 with permission from [32] p. 27. Copyright (2009) by the American Physical Society. 



production of W + n jets, a basic quantity to examine is the pr distribution for the 
softest observed jet. Fig. 2 shows the px distributions of the softest observed jet in 
W + ra-jet (n = 1,2,3,4) production at LO and NLO, respectively. For details on 
our analysis setup we refer to [36]. The predictions are normalized to the central NLO 
prediction in the lower panels. Comparing the pt distributions successively starting from 
W^-|-l-jet production, we observe the reduction in differential cross section of about a 
factor of as from one panel to the next; each observed jet leads to an additional power in 
the strong coupling. At the same time the lower panels in fig. 2 show an approximately 
linear increase in the LO scale variation bands, up to about ±50% for W + 4 jets. The 
scale variation of the NLO result, displayed in the lower panels of fig. 2, is strongly 
reduced to about ±10% for the present setup. 

In summary, in the above examples the advantages of NLO computations over 
the leading order appear through several quantitative improvements. Firstly, NLO 
predictions show a greatly reduced dependence on unphysical renormalization and 
factorization-scales as compared to leading order. The second improvement we pointed 
out concerns the shapes of distributions. Due to inclusion of radiation from an 
additional parton as well as a more truthful description of the scale dependence shapes 
of distributions are modeled better at NLO. 
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Figure 2. A comparison of thepT distributions of the softest observable jets in VF~+n- 
jet {n — 1,2,3,4) production, respectively. The setup describes the LHC proton-proton 
collisions at — 7 TeV as published in [36]. In the upper panels the NLO distribution 
is the solid (black) histogram and the LO predictions are shown as dashed (blue) lines. 
The lower panels show the LO distribution and LO and NLO scale-dependence bands 
normalized to the central NLO prediction. 



The above results validate our understanding of the W^-|-n-jet processes for typical 
Standard-Model cuts. It will be interesting, and necessary, to explore the size of 
corrections for observables and cuts used in new-physics searches. A related process 
that contributes an irreducible background to certain missing energy signals of new 
physics \s Z + 4-jet production. We expect that the current setup [36] will allow us to 
compute NLO corrections to 4-jet production, as well as to other complex processes, 
thereby providing an unprecedented level of theoretical precision for such backgrounds 
at the LHC. 

Parton-level NLO simulations of this kind are first principle predictions whose 
outcome directly reflect properties of the underlying theory. Although NLO 
computations are more challenging, in general they yield results with better reliability 
and agreement with measurements. 

2.2. Setup of Complete Computation 

The computation of differential distributions is the end product of combining many 
important ingredients pulled together in a Monte Carlo program; these include parton 
distribution functions and couplings, phase-space integration, matrix elements, analysis 
framework etc. Various tools are available to deal with complete NLO computations. 
One such tool is MCFM [108], which contains an extensive library of analytic matrix 
elements for NLO computations. Another approach (see [35] and references) uses tools 
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Figure 3. Representative diagrams of matrix elements for W + 4-jet production: 
(a) the eight-point loop amplitudes qg — evq' ggg, and (b) the nine-point tree-level 
amplitudes qg evq' gggg needed for the real contribution. The ev pair couples to 
the quarks via a W boson. 

including Helac [7] and CutTools [79] for a numerical approach [27]. Here we will 
describe another setup [30, 32, 36] based on on-shell and unitarity methods that was 
used for the computation oiW + n-jet production in section 2.1. 

In addition to LO components of Monte-Carlo programs, at NLO the computations 
rely on further similarly important ingredients. For the W^-l-n-jet production [30, 32, 36] 
the real-emission and dipole-subtraction terms [58], are provided by the SHERPA 
package [12]. SHERPA was used to perform phase-space integration. BlackHat was 
used to compute the real-emission tree amplitudes for Vr-|-4jets using on-shell recursion 
relations [47] , along with efficient analytic forms extracted from TV = 4 super- Yang-Mills 
theory [109]. 

In terms of scattering amplitudes we need the input of up to eight-point one- 
loop QCD amplitudes as well as up to nine-point tree-level QCD amplitudes; example 
Feynman diagrams are shown in fig. 3. The squared matrix elements are summed 
over all initial and final-state partons, parton helicities and color-orderings. For the 
present computation the W-hoson is decayed into an observable electron and a neutrino. 
Amplitudes of this kind can be obtained from QCD amplitudes; with the lepton pair 
replaced by a quark pair and the W-hosow exchange related to a gluon exchange. 
Appropriate dressing with coupling constant and propagator terms are needed. A recent 
analysis of high multiplicity tree amplitudes of this kind can be found in [109]. 

3. Structure of One-loop Matrix Elements 

The evaluation time of matrix elements is often dominating cross section computations, 
thus, emphasising the importance of efficient numerical algorithms. Beyond this, matrix 
elements are objects of fundamental theoretical interest; new physics effects observable 
at high energy colliders may originate in properties of matrix elements (see e.g. [110]). 

Matrix elements are functions of a large number of variables, which characterize 
particles, polarization states, color quantum numbers, and kinematics. To next-to- 
leading-order in the strong coupling the parton level cross sections for (A — 2) resolved 
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final-state objects, pp — > {N — 2), depend on squared born matrix elements, 

J2K''i{ki,hi,ai}\', n^N,N + l, (1) 

as well as the interference terms, 

J2^T*{{h:h,,a,})Al;"'^n{ki:h,,a,}) + c.c. , n^N, (2) 

where /cj, /ij, and Oj are respectively the momentum, helicity (±), and color index of 
the i-th external gluon or quark. The shorthand c.c. refers to the complex conjugate 
part that has to be included. The efficient management of parton and helicity sums is 
important. For simplicity, we will consider scattering amplitudes involving quarks and 
gluons only. Much of the methods can be carried over to more general particle spectra. 

As inspired by analytic approaches (see e.g. [Ill]) we disentangle degrees of freedom 
in order to arrive at a fine set of gauge invariant objects. To this end several structures 
are used: Color decomposition into color ordered sub amplitudes disentangles color 
information and kinematics. Use of spinor helicity notation aligns notation of kinematics 
and polarization vectors. Spinor variables, in addition, lead to a natural way to work in 
complex momentum space. This in turn allows to exploit analyticity properties of the 
basic color ordered scattering amplitudes. Use of a standard basis of integral functions 
will allow a further fine split-up of the loop amplitude into integral functions and their 
integral coefficients. 

Several features motivate us to disentangle matrix elements into a fine set of 
primitive amplitudes. First of all, do these have cleaner physical and analytic properties 
than the full matrix elements, as will be discussed in detail below. This can be exploited 
for the construction of computational algorithms, as will be discussed in detail below. 
For numerical approaches, a detailed understanding of physical properties (e.g. IR/UV 
pole structure of primitive amplitudes, or, tensor rank of integrals) allows one to monitor 
precision and stability of the computation. Fiuthermore. caching systems built on 
primitive objects (here, amplitudes) lead to important efliciently gains through reduction 
of redundant computations. Finally, during numerical phase-space integration, one can 
introduce importance sampling, computing computationally expensive but numerically 
small parts on a subset of phase space points. 

3.1. Color Decomposition 

At tree-level, amplitudes for SU{Nc) gauge theory with n external gluons can be 
decomposed into color-ordered partial amplitudes, multiplied by an associated color 
trace (see e.g. [74]). Summing over all non-cyclic permutations gives the full amplitude, 
AT{{h,h,,a,})= Tr(TMi)...rM.))^;^-(^^.a)^_ _^^V.)) ^^3) 

the coupling is set to one, and Sn/Zn is the set of non-cyclic permutations of {1, ... ,n}. 
The are the set of hermitian generators of the SU (Nc) color group. The coefficients of 
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the color structures Tr(T"i • • • T"") define the tree-level color-ordered partial amplitudes, 
A^:-(l,2,3,---n). 

One of the important features of this set of amplitudes is that it forms a 
closed set under collinear, soft and multi-particle factorization. They have manifest 
transformation properties under parity transformation and reversal of the ordering of 
external legs. Similarly, amplitudes with fermions can be decomposed into color-ordered 
sub- amplitudes [112]. 

For one-loop amplitudes, one may perform a similar color decomposition to the one 
at tree- level in eq. (3). In this case, there are up to two traces over color matrices [76], 

\n/2\ + l 

^^-i-P({A;,,/i,,a,})= E Gr„;e (a) (4) 

where \_x\ is the largest integer less than or equal to x. The leading color-structure 
factor Gr„.i(l) = Tr (T"^ ■ ■ ■ T"" ) is just N^. times the tree color factor, and the 
sub-leading color structures arc given by the double trace expressions, Gr„;c(l) = 
Tr (T"i • • - T"-"-^) Tr [T""" ■ ■ -T""). Sn is the set of all permutations of n objects, and 
Sn-c is the subset leaving Gr„;c invariant. 

The leading partial amplitudes A^-i take the form (see e.g. [113]), 

<1°°P(1, 2, 3, • • • , n) = yl^(l, 2, 3, • • • , n) + ^A^(l, 2, 3, • • • , n) , (5) 

with Af^ and being color-ordered sub-amplitudes, or primitive amplitudes. While Af^ 
is fixed to have only gluons propagating in the loop, is restricted to have a Weyl 
fermion propagating in the loop. The external gluons are color-ordered; the diagrams 
contributing to the loop amplitudes can be generated from color-ordered Feynman rules, 
see e.g. ref. [113]. 

The coefficients of the sub-leading color structures; the sub-leading partial 
amplitudes, can be expressed in terms of sums [77, 78, 76] of the primitive amplitudes, 
Af^ where different ordering of the external states appear. Beyond the fact that such a 
decomposition exists, we will not need details of its form here. 

The primitive amplitudes, 

Ag(l,2,3,---,n) and ^1(1, 2, 3, • • • , n) , (6) 

form a generating set of amplitudes, such that given these amplitudes, the full one-loop 
matrix elements can be computed. For fundamental fermions a similar split-up of partial 
amplitudes is typically more complicated [78, 24]. In addition to the ordering of the 
external leg the routing around the loop (left- or right-turner) of each of the fermion 
lines has to be specified. See figure 4 for examples of primitive amplitudes including 
also external fermion lines. 

Primitive amplitudes, have a more transparent analytic structure than full matrix 
elements, because their external colored legs follow a fixed cyclic ordering. In particular, 
properties of factorization and branch cut singularities simplify as can be summarized 
by the following: 
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Figure 4. Representative Feynman diagrams for primitive amplitudes with distinct 
states propagating in the loop: (a) shows a gluon loop of the primitive amplitude Af, 
(b) a fermion loop of and (c) a mixed fermion/gluon loop of ^4^. For the mixed 
amplitudes (c) we keep track of the routing of the fermion line around the loop; 'L i?' 
indicates that the first fermion lines turns left and the second right when entering the 
loop. 



(i) Only factorization poles and branch cut singularities in adjacent legs appear. 

(ii) Primitive amplitudes and color-ordered tree amplitudes form a closed set under 
factorization and unitarity cuts. 

With the notation specification of 'adjacent legs' we refer to the fact that factorization 
poles and unitarity cuts appear only in a specified subset of kinematic invariants 
Sij^...if., with an ordering of momenta identical to the one of external gluons. Closure 
under factorization means that color ordered amplitudes factorize onto color ordered 
amplitudes. In particular, for unitarity approaches as well as on-shell recursions, 
primitive loop amplitudes can be constructed from color-ordered tree amplitudes alone. 

For a more complete description of color decomposition we refer the reader to 
previous reviews [78, 24, 74]. 

3.2. Structure of Loop Amplitude 

Any one-loop amplitude can be written as a sum of terms containing branch cuts in 
kinematic invariants, C„, and a rational remainder 

An = Cn + Rn ■ (7) 

The cut-containing part C„ can in turn be written as a sum over a basis of scalar master 
integrals [114, 115], 

«l<i2<*3<«4 «l<i2<«3 «1<«2 

The scalar integrals 12^3^4 - bubbles, triangles, and boxes - are known functions [116]. 
We omitted tadpole functions, which in dimensional regularization vanish for massless 
particles circulating in the loop. The explicit dimension dependence is contained in the 
integral functions with their coefficient functions strictly four dimensional. Feynman 
diagrams of the integral functions are shown in fig. 5. 
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Figure 5. Feynman diagram representation of basis of scalar integrals appearing in 
the eq. (8): (a) box diagram associated to the 4-point integral, (b) triangle diagram 
associated to the 3-point integral, and (c) bubble diagram associated to the 2-point 
integral. Each corner may have one or more external lines attached to it. 



The computation of a one-loop amplitude amounts to determining the rational 
coefficient functions c° and 6°| in addition to the remainder Following the 

spinor-helicity method [73, 81], we can then re-express all external momenta in terms of 
spinors. The coefficients of these integrals, c°, and 6°, as well as the remainder i?^, 
are then all rational functions of appropriate spinor and momentum variables. 

For the analysis of one-loop amplitudes it is often useful to have two distinct forms 
of the integrals in mind. One can think of the integral functions as logarithms and 
polylogarithms of kinematic invariants. As examples we give explicitly a bubble integral 
function, 

X2(,i2, fi^) = icr{\+{2 + ln(3^) } + 0{t) , (9) 

exposing discontinuities in kinematic invariants through branch cut singularities of 
the underlying logarithmic functions. For the kinematic invariants Sij the standard 
definitions, Sij = {ki + kjY, are used. The constant cr is defined by, cr = 
r(l -t- e)r^(l — e)/((47r)^~'^r(l — 2e)) . The scale fi corresponds to the renormalization 
and factorization scales which, for convenience, are set to be equal here; fi = fiR = fip- 
The integrals are also viewed as Feynman amplitudes of a scalar field theory, 

where sums of two momenta are denoted by the shorthand, Kij = ki + kj. 
3.3. The Loop Integrand 

For explicit computations it is useful to consider loop amplitudes before integration, 
that is to find a universal parametrization [27, 29] of the loop integrand. In addition 
to the scalar Feynman diagrams fig. 5, implied from eq. (8), tensorial numerator terms 
have to included. These tensorial terms describe as well the angular distribution of the 
virtual particles, which averages out upon integration. The explicit relation between 
numerator tensors and angular variables will be topic of section 4.3. 

I As a shorthand we often specify the collections of indices {zi, 12, . . .} by a single one; e.g. instead 
of ■ 
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A generic form of the loop integrand is given by [27, 29] (see also [28]), 




(11) 



where d stands for a given discrete dimension and we restricted external momenta and 
polarizations to be four-dimensional. If we allow d-dimensional polarization vectors and 
external momenta, higher polygons have to be added in a natural way. The pentagon 
terms e{€) should be dropped when working in strictly four dimensions, d = 4. In the 
above expression, propagators are denoted by i/Dj — ijii — KjY\ for simplicity, of 
notation we restrict the discussion on massless internal states. Furthermore, we omitted 
the single propagator terms a(£) which drop out of the final results in the absence 
of massive states. When used with the explicit argument i for the loop momentum, 
74^~^°°P''^(£) denotes the integrand as opposed to the loop amplitude 

The numerators e(£),(i(£), c(£), h{t) and a(£) are sums of tensors contracted with 
the loop momentum l^. The tensor-rank is bounded by power-counting. We will refer 
to these tensors contracted with loop-momentum I somewhat imprecisely as numerator 
tensors. These numerator tensors can be expanded in terms of a basis of tensors in 
momentum space multiplied by scalar loop-momentum independent coefficients. The 
scalar coefficients then characterize the loop amplitude. See below in section 3.3.3 for 
an explicit representation in terms of a basis of tensors and scalar coefficients. 

Integrand parametrizations (11) are common in unitarity approaches; for a 
discussion in the context of multi-loop computations sec e.g. rcfs. [71]. A particularly 
useful parametrization of the one-loop integrand has been given in [27, 28, 29], as will 
be discussed further in the following section. 

3.3.1. Loop Integration. With an appropriate representation of the loop integrand the 
loop integrations can be performed trivially. This is achieved by writing the integrand 
numerators as a direct sum [27] of terms that integrate to zero and non- vanishing scalar 
terms. The form of the integrand in eq. (11) then directly relates to the from eq. (8). 
When we evaluate the known analytic expressions for the basis of integrals, we thus 
obtain an exact numerical algorithm to go from an off-shell integrand to the integrated 
loop amplitude. An approach of this kind was used in refs. [27, 21, 28, 29]. We will 
motivate a canonical form for the loop integrand in this section leaving a more complete 
discussions to the next. Such a choice can be viewed as an implicit integral reduction 
procedure. 

To have an example in mind, consider the box numerators in the form. 



Here the vector is understood to be orthogonal to the external momenta of the 
box function. (For an explicit definition and properties of n]^ see later in eq. (17).) 




rO 



'1 



1 



(12) 
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The coefficients and dj ■ ■ ■ are the free parameters of the ansatz that have to 

be determined. The coefficients we eventually need to compute are the coefficients 
of the scalar term d^_^^,^,j^,^j^^ which correspond to the scalar basis integral coefficients 
in eq. (8). The tensor coefficient d}_^i^i^i^, although necessary at intermediate steps 
of the computation drops out after integration, 

(^^■^) (13) 

That is, after integration the numerator loop momentum gets replaces by a linear 
combination of the external momenta Ki which are orthogonal to n^; giving a vanishing 
tensor integral. 

It is instructive to consider another form of the tensor numerator, including a term 
with an inverse propagator, 

di,i,isM = <,,,3,, + <,,,3,,(n^ • e) + <,,,3,,A,(^) . (14) 

Clearly the inverse propagator may be cancelled against one of the box propagators 
turning d^^^^^^^^ into the coefficient of a scalar triangle function. If we were to treat this 
term at the box level, we obtain tensorial contributions that have to be transcribed into 
scalar integral functions with some care. Assuming for the moment some prior knowledge 
about generalized unitarity cuts we can make some further observations. That is, this 
term actually vanishes on the quadruple cut, leaving the coefficient (^f^j^irU unspecified 
at first. It may then be fixed using the triple-cut equations, although only the sum of 
the scalar triangle and d1_^^^j^^^^ may be fixed. The numerical unitarity relations eq. (51) 
are then not triangular. Then box, triangle and bubble coefficients cannot be solved for 
consecutively. This is of course obvious in the present example, given that we rewrote a 
scalar triangle coefficient in box-form. This observation emphasises the need for a good 
basis of numerator tensors. 

A less trivial deformation of the numerator tensors would be to mix in a propagator 
term with the linear tensor, 

d^,^,^,M = <,,,3,, + <,,,3,, (Ai(^) + {u^ ' t}) . (15) 

Here again a scalar triangle contribution is pulled back into the box integral. This time 
the triangular nature of the cut-equations (see below in eq. (51) for the explicit form 
of the cut-equations) stays intact and no redundancy is introduced into the numerator 
tensors. However, one has to pay attention not to drop the coefficient when integrating 
the tensor box integrals. We can read off the box coefficient directly, d'^^^^^^^^, and in 
addition, we have to include d]^^^^^^^ to the related triangle coefficient. 

Terms may be moved around between integral functions in this way, effectively 
introducing a change of basis of integral functions. As a non-trivial application, a cut 
completion, that is a subtraction of gram-determinant poles, may be achievable in this 
way. The form of the numerator tensors has important implication; it allows to keep the 
unitarity relations triangular, and, keeps the integration of the loop integrand simple. 
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3.3.2. Numerator Tensors. In the numerical unitarity approach one is naturally lead 
to obtain equations for the integrand expression eq. (11). It is then convenient to use 
an explicit form of the numerators in terms of a basis of tensors. Computing a loop 
amplitude then amounts to determining the free tensor coefficients. 

There arc several natural requirements [27] for a good basis of numerator tensors. 
The first requirement is, that numerator tensors should of course be general enough 
to parametrize the loop-integrand we are interested in. Typically one uses all tensors 
up to a given rank, as determined by power-counting. Furthermore, optimally one 
would like to use a minimal set of tensors. A last requirement is then, that it should 
be easy to relate the integrand basis back to the integral representation in eq. (8). 
It turns out that an optimal tensor basis can be found, which satisfies all the above 
requirements [27, 28, 29]. 

For numerator tensors in strictly d dimensions the tensor basis looks particularly 
simple. (We will discuss the D-dimensional generalizations below in section 4.8.) In fact, 
the result will be a basis of tensors, called spurious numerators in the literature [27], 
which integrate to zero, 

{27rr-^^e{i-K,y...{i-K,...-K^.,r^ ^^^^ 

where n^^...^^ stands for a representative basis tensor. Upon integration, the loop- 
momentum dependent numerators in eq. (11) may thus be dropped and the remaining 
scalar (rank-zero) terms are directly identified with the integral coefficients in (8). 

In order to obtain this basis of tensors it is convenient to introduce the Neerven- 
Vermaseren basis [115] for vectors in momentum space; a distinct basis for each of the 
integral functions. Each integral defines a distinguished set of momenta; the momenta 

in eq. (16). Momentum space is decomposed into the direct sum of two subspaces; 
the physical space parametrized by and its complement, spanned by the vectors 
n'^" [28], 

Ki for i e {1, • • • ,n - 1} , for i e {1, • • • , d - (n - 1)} , 

■ = 5'^ , n'-Kj = 0, 
v' for i G {1, ■ ■ ■ , n - 1} , Ki ■ = Sj , n' ■ = , (17) 

where we assumed n < d. For the complementary case n > d momentum space is 
parametrized solely in terms of a linearly independent set of vectors K^. 

The vectors v^'^ are dual to the external momenta § and are part of the physical 
space, defined by the external momenta K^^. The vectors n*'' are an orthonormal basis in 
transverse space and are orthogonal to the physical space. Depending on the signature 
that the transverse space inherits from momentum space, the vectors n*^ have to be 
chosen purely real or imaginary. A further useful quantity is the metric g'^ of the 

§ We may obtain the dual vectors using = Ki ■ Kj and it's inverse {k~^)''^ , such that = 
^ ■(k;"^)*'' /^^. An exphcit form of the dual basis can be found in [28]. 
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transverse space, 

{gi_)^''^' = n'^^n'^' = ^^^^ - ^ v'^'K^^ , (18) 

i=l i=l 

which is naturally related to the metric of momentum space. 

A generic numerator tensors can be expressed as tensor products of the vectors (17). 
A basis of tensors is thus given by, 

n,.-,, =Kl^...Kii^r^,...rvlz, l + m = k. (19) 

In fact, the set of all numerator tensors needed is given by the symmetric traceless 
tensors in the transverse space, 

^™ (20) 

where the curl- brackets denote the operation of symmetrization and subtraction of traces 
within transverse space. By trace we mean the contraction, 

n,.,,-,^9l''' = , (21) 

of Lorentz indices with the metric tensor in transverse space, g'^. For symmetric tensors 
it is sufficient to single out two indices for contraction. 

Again the case of a vanishing transverse space n > d is special. For this we are left 
only with the rank-zero scalar numerators. In fact, for this case a further simplification 
appears; the scalar n-point integral can be written in terms of lower-point integrals. 
To show this one uses identities implied by inserting vanishing Gram determinants, 
A(i+i(iri, • • • Kfi,i) — 0, into n-gon integral. Repeated reasoning along these fines leads 
to reduce n-gon integrals to d-gons or lower. We refer to a recent discussion on this 
in [117] for further details. 

Examples of symmetric traceless tensors are then, n^^rh) = {nj^nl + n'j^nl)/2, which 
is traceless due to the orthogonality of the vectors n*. A further example is the tensor, 
njt^nD = {njj^nl — {g±)fj,u/{d — n + 1)), for which the trace was explicitly subtracted. 

The form (20) of the spurious terms can be understood in the following way. To 
start with, it turns out that tensors (19) with components pointing along the physical 
space, i.e. I > 0, are redundant. For the simplest case — (with I — 1) the 
contraction of Ki with the loop momentum £, 

i.K, = ^[e-{i-K,)' + Kl] , (22) 

gives rise to inverse propagators £^ and {£'^ — KiY, and a scalar term (Kf). Although 
we started with a rank-one tensor integrand, after the inverse propagators is cancelled, 
we obtain lower-point integrals and a scalar integral, 

Ki-i 1' 



+ 



{£ - Xi)2L>2 ■■■Dn £^D2 ■■■Dn 

£^{e-K,)W2---D^ 



+ 



(23) 
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The tensor integrand we started with is thus redundant, as it can be expresses solely 
in terms of lower-rank and lower-point terms. A similar reasoning, applied recursively, 
shows the redundancy of tensors with multiple fCj-components and enforces / = in the 
notation of eq. (19). Tensors including physical directions K-i are thus redundant; they 
either lead to linear combination of lower-point integrals or are expressed as tensors of 
lower rank. We thus do not need to consider these tensors further, once we account for 
this. 

What remains to be considered are tensors with components purely in the transverse 
space. Of these only a subset of tensors is linearly independent. In particular, a trace- 
containing term in the transverse space is related to a trace in physical space in addition 
to a metric tensor. Thus, a trace-containing term yields inverse propagators and loop- 
momentum independent terms when contracted with loop momentum, 

d—n+l n—1 

^^^'^ = - ■ ^)(^^ ■ ^) ' (24) 

i=l 1=1 

where we used equation (18). The contractions of {Ki-i) can be transcribed into inverse 
propagators and terms independent of loop momentum as in eq. (22). Traces thus lead 
to lower-point or lower-rank tensors and are thus linearly dependent. The only choice 
left are traceless, symmetric tensors in the transverse space. 

One might still worry that additional hidden relations can be found to relate 
integrals with distinct tensor numerators of the from (20). No further relations can 
in fact be found. The independence of the tensor of the same propagator structures 
can be argued using an explicit on-shell loop-momentum parametrization in eq. (52) 
as we will discuss further in section 4.3. The independence of tensor integrals (20) 
with distinct propagator structures is due to their differing factorization properties; 
e.g. triangle integrals with numerator tensors (20) cannot mimic the quadruple cut 
divergence of four-point integrals. 

Tensor integrals with symmetric traceless numerator tensors integrate to zero. Due 
to Lorcntz- and parity-invariance, see e.g. [27, 118], a generic tensor integral is written 
in terms of productions of the vectors and metric tensors g^v, 

^ J (27r)4-2^ P{i - K,y ...(£- ... - ir„_i)2 - 

= kI^' . . . . . . g^'^-^^^'^^fisij) + ■■■ . (25) 

The integrals simply depend on no other vectors and e-tensors are excluded by parity 
invariance. Upon contraction with a symmetric traceless tensor, n^^...^^ from (20), the 
left hand side of (25) turns into a tensor integral (16) while the right hand is easily seen 
to vanishes; n^^...^^ is traceless and the vectors n*^ being orthogonal to K^. We thus 
verified that symmetric traceless numerator tensors (20) lead to vanishing integrals (16). 

In summary, the symmetric traceless tensors (20) fulfill all criteria of an optimal 
basis as discussed in the beginning of this section. 
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3.3.3. Tensor Basis. An explicit form of the numerator tensors in eq. (11) in terms of 
the vectors (17) was given in [28]. For box coefficients we have, 

for triangles, 

Cni2i3(^) = (^111213 + ^\ii2V.i^^ + ^\i2kM + ^\i2i?X^^ ^1 ~ ^2 ^2) 

+ ^1 h {cli,i, + cl,^i,h + cl,^i,h) , (27) 

and bubbles, 

+ hU,{hh-hh) + hl,^{hh-hh) 

+ hl,.trt2 + hl,^Uh + hl^^t2h, (28) 

respectively. Here we introduced tj = (n* ■ I). The vectors Uj differ between the 
three equations and arc defined for each associated propagator structure individually as 
defined in eq. (17). The coefficients we wish to compute arc the and 6?,,- 

J-V/ ^ ^l^z''3^4: '1^2''3 ^1''2 

terms which correspond to the scalar basis integral coefficients. The tensorial expressions 
vanish upon integration but have to be kept at intermediate steps of the computation. 

The above representation is not unique; not only may one chose a different basis for 
the transverse space and thus different basis vectors . One may also alter the tensor 
basis used in this parametrization. For example, the above expression uses terms of the 
form (ti ti 12) which are not tracclcss. An alternative, traceless symmetric representation 
would be instead (titit2 — ^2^2^2/3). The difference of the two approaches amounts to 
moving loop-momentum tensors between bubble- and triangles-functions. Given that 
both forms integrate to zero, either of the above choices leads directly to the same final 
result. (For a related discussion see also section 3.3.1.) 



3.4. Spinor Helicity 

Spinor variables give a unified way to express polarization vectors of gluons, fermion 
helicity states and kinematics of a scattering process. Furthermore, spinor variables lead 
to a natural way to work with complex momenta. Complex-valued on-shell momenta 
are important in order to fully exploit analyticity properties of amplitudes. We will see 
examples of this in computations of integral coefficients, sections 4.4 and 4.7, and, later 
in section 5, when we consider on-shell recursions. 

We follow the standard spinor helicity notation and conventions as in refs. [73, 74] . 
As a shorthand notation for the two-component (Weyl) spinors we use, 

(Xi)^ = [u+(ki)]^ , {Xi)^=[u4ki)U, kf^O. (29) 

Lorentz-covariant spinor products of left- and right-handed Weyl spinors can be defined 
using the antisymmetric tensors e'^^ and e°^, 

{j I) = e^'\\M^i)p, [j I] = e"^{XM~^i)3 ■ (30) 
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These products are antisymmetric, {j I) = — (Ij), [j I] = — [Ij]- 

One can reconstruct the momenta from the spinors, using u{k)u{k) — ^, 

^i{'^n)aa — {l^i)aa — {\)a{\)a ■ (31) 

Equation (31) shows that a massless momentum vector, written as a bi-spinor, is simply 
the product of a left-handed spinor with a right-handed one. In order to specify on-shell 
momenta we will often use the abbreviation, 

£-Ai(^A,), (32) 

where spinorial indices are suppressed and the index contractions are indicated by the 
parenthesis; {^\j)a = ^aa^'j- Spinor products of the above momentum i are then 
given by, {ni) = {ni) and [in] = {j\K\n]. 

The usual momentum dot products can be constructed from the spinor products 
using the relation, 

{lj)[jl] = 2kyki = sji. (33) 
We will also use the notation, 

m 

{a\Kn...iJb] ^ ^ ^n-im ^ Kl,„i^ . (34) 

k=l 

A further important class of quantities are Gram determinants An{Ki, ■ ■ ■ , Kn) 
defined by, 

A„(Xi, ■■■,Kn)= det(2 Kj ■ Ki) . (35) 

Gram determinants appear naturally in unitarity cuts; when solving for on-shell 
momenta negative powers of Gram determinants appear. These then enter the 
computation of integral coefficients. 

Gram determinants can be associated to the linearly independent momenta of 
integral-functions. The respective integral coefficients typically have inverse powers 
of these Gram determinants in addition to the ones inherited from reduction of higher- 
point tensor integrals. For loop-level on-shell recursions, section 5.3.2, we will see that 
Gram determinants play an important role. 

3.4- 1- Basic Tree Amplitudes. When using the spinor-helicity formalism tree-level 
scattering amplitudes simplify significantly. Further simplifications arises in part from 
symmetry properties of tree amplitudes [50, 51, 52] that are present in QCD-like theories 
(see also [53]). Numerical implementations of on-shell recursions may recurse all the way 
to three-point vertices. More efficiently, they can easily be combined with a library of 
compact analytic trees. The recursion is stopped and analytic expressions are used 
whenever available, leading to an efficient numerical algorithm. 

One example of simplifications due to the use of spinor helicity variables are the 
infinite set of vanishing tree-level gluon amplitudes, 

Ar(1^2+,...,n+) = 0, n>3, (36) 
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with all helicities identical, or all but one identical. Parity may of course be used to 
simultaneously reverse all helicities. 

The infinite set of Parke- Taylor amplitudes [119, 120] is another striking example 
for which the use of spinor helicity formalism yields a particularly simple form, 

with two negative helicities and the rest positive. The only gluons with negative helicity 
are in positions 1 and j. Helicities are assigned to particles with the convention that 
they are outgoing. The parity conjugate amplitudes may be obtained by exchanging 
the left- and right-handed spinor products in the amplitude, {j I) ^ [I j]. 

Furthermore, implicit supersymmetry properties [121] allow to relate fermion, gluon 
and scalar amplitudes of differing spins. For example, in order to replace the gluons 1 
and n by scalar states, 

^tree(^- 2+ J- ^.^^+) j^;^') ri > 3 , (38) 

^ ^' ' '-^ ' ' (1^)2 (12) (2 3)- --(nl) - ^ \ J 

we simply multiply the pure gluon amplitude with an overall factor. Such relations can 
be used to speed up the sums (1) over final-state particles when using trees with manifest 
supersymmetry properties . (See e.g. [49, 122] for trees with manifest supersymmetry 
properties.) 

3.4-2. On-shell Momenta. For real momenta, Aj and Aj are complex conjugates of each 
other up to a sign depending on the sign of the energy component. For the degenerate 
but important case of three-point kinematics, 

h + k2 + h = ^, with A;f = , (39) 

only the trivial solution ki ~ kj can be found for real momenta. For these real solutions 
all spinor products vanish. 

However, for complex momenta, it is possible to choose all three left-handed spinors 
to be proportional, Ai = C1A3, A2 = C2A3, while the right-handed spinors are not 
proportional, but obey the relation, ciAi -I-C2A2 + A3 = 0, which follows from momentum 
conservation, ki -\- k2 -\- k^ — 0. Then, 

(zj)7^0, but [zj] = 0. (40) 

A second branch of solutions to the on-shell conditions can be found as the conjugate 
set of momenta, Aj ^ Aj. 

Such degenerate kinematics are important for unitarity cuts associated to integral 
functions with massless corners. An explicit computation will be discussed in section 4.4 
and section 4.7. For such cases three-point tree amplitudes, 

(12^^ 

^(1-, 2-, 3+) = I , ,\ ; , , , (41) 
^ ^ ' ' ^ (12)(2 3)(31)' ^ ' 
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have to be evaluated on solutions to the on-shell conditions (39). They are non-trivial 
on one set of complex solutions and vanish as 0^/0^ on the other. The above non- 
trivial solutions involving complex momenta are then necessary in order to exploit 
generalized unitarity cuts. The general form of this type of on-shell conditions is 
discussed below eq. (55). 



3.5. Supersymmetric Decomposition 

The supersymmetric decomposition of the amplitudes is particularly useful when 
considering rational terms of scattering amplitudes. In particular, the rational parts 
of amplitudes with gluon and fcrmion degrees of freedom circulating in the loop can be 
related to often easier to compute scalar ones. 

From power-counting arguments we known a priori that the supersymmetric 
amphtudes A^^'^ and A!^^^ are cut constructible in four dimensions and are free 
of rational terms [22]. The M — ^ multiplet and the M — 1 chiral matter multiplet 
are built from a particular combination of gluon, fermion and scalar degrees of freedom. 
For the case of external gluons, the couphngs of the matter particles resembles the one 
of QCD, leading to the relations, 

A^=' = Ai-r 4Ai + 3AI , 

^Ar=l chiral = ^/ + (42) 

between supersymmetric amplitudes (Ihs) and basic field theory amplitudes (rhs). The 
superscripts, g, f and s, indicate the states circulating in the loop, and stand for gluon, 
Weyl fermion and a complex scalar, respectively. Although the above relations are for 
adjoint fermions in the loop, they can be directly related to massless fundamental quark 
loops [123, 22]. 

Inverting the above relations (42) one obtains the amphtudes for QCD via 

A9 — 4 A/'=4 _ ^ A/'=l chiral , as 

Ai= A„^=^'^^-^^ - a:. (43) 
This then implies that the rational terms within A^ and A-^ equal the ones from , 

I rational I rational : rational I rational ■ (^'^) 

With this decomposition we can then compute the cut containing pieces in strictly four 
dimensions taking into account the full QCD spectrum in the loop. At the same time 
one may compute the rational part of the QCD amplitudes purely from amplitudes with 
a complex scalar in the loop. 

When computing the rational terms using the D- dimensional unitarity ap- 
proach [23], virtual scalars are much more straightforward to deal with as opposed 
to gluons or fermions. While the kinematics of internal scalars has to be considered 
beyond four dimensions, we do not have to worry about D-dimensional extension of po- 
larization states of gluons and fermions. Computations are then very similar to having 
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a massive scalar [23, 124, 72] in the loop, however, where the mass is related to the 
[D — 4)-dimensional momentum. 

Relations of this kind are rather generic and can be found for internal fermions and 
mixed fermion and gluon amplitudes. For analytic computations this observations were 
used for example in refs. [78, 24]. 

4. The Unitarity Method 

The modern unitarity method [22] in addition to generalized unitarity [24, 25] are 
the foundation of powerful approaches for loop computations with phenomenological 
interest. Many recent generalizations [27, 28, 21, 29], in particular, with a numerical 
application in mind, have helped to established a standard unitarity algorithm. 
These numerical unitarity methods were first applied to studies of hadron collider 
physics in [30, 31, 32, 69, 34], and are by now used by many groups [37]. Beyond 
this, various other implementations of numerical unitarity approaches have been 
reported [40, 39, 41, 38, 42]. 

Below we will focus on key developments of the unitarity method with emphasise 
on numerical aspects. We will follow aspects of the approaches outlined in [28, 21, 29]. 
For a discussion of analytic unitarity methods we refer to the chapter in the review by 
Britto [43] and references therein. A more detailed account of the modern unitarity 
approach as well as its application for multi-loop computations may be found in the 
chapters of this review by Bern and Huang [125] and by Carrasco and Johansson [71]. 

4.1. Unitarity Relations 

In terms of the non-forward part T of the S-matrix, S = {1 +iT), unitarity conditions 
= 1, imply the nonlinear equations, 

-i{T - T^) = T^T. (45) 

When combined with analyticity properties, as present in field theory, the unitarity 
condition (45) relates branch cut discontinuities of scattering amplitudes, to integrals of 
products of scattering amplitudes. (See e.g. refs. [126] for an early account of unitarity 
and analyticity.) At one- loop order the unitarity relations may be written as, 

Disc(s,) = E / dmj2) A':ri-£2ji)A':^%-hj2) , (46) 

states 

where the state sum is over all intermediate physical states in the theory. The phase- 
space integral J ^$(^1,^2), is defined over integration contours with the intermediate 
momenta £1 and £2 on-shell; £f — mf. The notation, Disc(s), stands for the branch 
cut discontinuity in the complexified variable s. E.g. for a logarithm we have 
Disc(s) ln(s//i^) = 27ri such that the operator Disc(s) picks out the coefficient of the 
logarithm. Importantly, the nonlinear unitarity relation links on-shell amplitudes of 
different loop-order in perturbation theory. For simplicity, we restrict our discussion to 
color ordered amplitudes. 



Susy Theories and QCD: Numerical Approaches. 



24 



For field-theory amplitudes Cutkosky [127] generalized eq. (46) further, providing 
a prescription to directly compute more generic discontinuities [128]. An early version 
of generalized unitarity for generic field theories was demonstrated in [24, 25], including 
massless states and an arbitrary number of external particles. Specialized to one-loop, 
the discontinuities are given by phase-space integrals of multiple on-shell scattering 
amplitudes, 

Disc(s,„-..s,J^^-^-P= J] I ci$<-<----<-. (47) 

states 

As above, the state sums run over all physical states in the theory. Intermediate 
momenta are integrated over appropriate contours of their simultaneous on-shell phase- 
space. 

For one-loop computations generalized cuts include single-, double-, triple-, 
quadruple- and penta-cuts. The generalized unitarity cuts are important for several 
reasons: first of all, the additional unitarity relations give further equations to 

characterize loop amplitudes. Furthermore, the additional on-shell conditions for the 
intermediate momenta restrict the phase-space integral further. The cuts that are the 
easiest to evaluate are the maximal cuts. For these the loop-momenta are frozen to 
specific values and the phase-space integral degenerates to only a sum over discrete 
solutions of the on-shell conditions. In four dimensions these correspond to quadruple 
cuts [25] . Finally, generalized unitarity allows to order cuts hierarchically; from maximal 
cuts to next to maximal cuts etc. This hierarchy allows a systematic approach to 
inverting unitarity relations. (See below in section 4.2.) 

Applied to Feynman diagrams, the generalized unitarity relations can be formulated 
as diagrammatic cutting rules [127]. In this setup cutting replaces propagators with on- 
shell conditions, 

ran; - ^ m 

and yields directly the values of the associated discontinuities. (The subscript 'p' in 
5p{l'^—m?) indicates the common restriction to a specific branch of on-shell momenta.) In 
this way cutting rules (48) relate unitarity cuts (47) to universal factorization properties 
of the loop integrands. More explicitly, under the loop integral propagators are replaced 
by on-shell conditions, such that conditions for the factorization limits of the loop 
integrands are obtained, 

hm (A,---A.Ar'°°^(^)) = 

(<n4, . . . , -4) X • • • X <^n4-i, . . . , -4)) . (49) 

states 

where A^~^°°''(£) stands for the loop integrand of the amplitude A^"^°°p. The state 
sums run over the full spectrum of the theory. The momentum solves the on- 

shell conditions Di^ = . . . = Di^ = 0, restricting the loop momentum to the on-shell 
phase-spaces. All intermediate momenta £i are on-shell and related to i by momentum 
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conservation. The importance of the use of complex momenta for solving the on-shcll 
conditions for theories with massless states was pointed out in ref. [25] where a algebraic 
equations for quadruple cuts were obtained similar to the form of eqn. (49). 

4-2. Inverting The Unitarity Relations 

We specialize our discussion below on computations in strictly four dimensions. 
Analogous ideas generalize to higher dimensions; see section 4.8. In principle, dispersion 
integrals [126] should allow us to assemble the amplitude in terms of its branch 
cut structure, however, for practical applications a more powerful approach can be 
used [77, 24]. The unitarity relations may be most directly implemented by relying 
on the decomposition of loop amplitudes into a basis of loop- integral functions (8). 
Matching the unitarity cuts with the cuts of basis integrals provides an effective means 
for obtaining expressions for the integral coefficients di, cj and bk, 

Disc A^"^°°P ^^^^ ^l+Yl ^^^^ ^3 ^^^^ ^2 ■ 

i j k 

The unitarity cuts of the ansatz are then compared in all channels to the cuts of the 
amplitudes. 

For a numerical algorithm one aims to further simplify eq. (50) by removing the 
phase-space integrations. This can be achieved by introducing a parametrization [27] 
of the loop integrand (11). In addition to scalar integrals, tensor integrals have to be 
considered. The latter parametrize the angular distribution of virtual states circulating 
in the loop, as discussed above in section 3.3. 

The computation of loop amplitudes then amounts to determining the coefficients of 
a basis of scalar and tensor integrands (26, 27, 28) and subsequently performing the loop 
integrations. Coefficients in the ansatz of (11) may be determined through comparison 
to factorization limits of loop amplitudes (49). The unitarity method applied at the 
integrand level [26, 28] leads to a set of equations for on-shell loop momenta, 

states 

states j^ii, 12,13 



states j^ii,i2 ^ j,k^ii,i2 ^ ^ 

where the individual equations have to be evaluated on respective on-shell momenta; 
i — >■ ini2i3uiz or ^ — >■ ii-^i2i3i4^, etc. The momentum dependence of the individual 
tree amplitudes is indicated with the appropriate shift of the momentum dependence 
implicitly assumed. 

Typically, the on-shell conditions do not constrain the loop momenta completely, 
but lead to a variety of solutions. Enforcing the unitarity relations (51) on the 
variety of on-shell solutions gives an infinite set of equations. The set of tensorial 
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structures in addition to the scalar integral coefficients can be determined due to this 
degeneracy [27, 26, 28, 21]. 

4-3. Angular Dependence of Numerator Tensors. 

It turns out, that the unitarity relations (51) allow us to obtain the coefficients of all 
scalar and tensor numerator terms, and, thus, allow to reconstruct the off-shell loop 
integrand. This is despite the fact that the unitarity relations are defined only for on- 
shell intermediate momenta. We will discuss a key aspect of this off-shell power of the 
on-shell unitarity equations in this section. 

An important observation is, that the on-shell conditions from cutting propagators 
fix the values of the loop momenta in the physical space, but introduce almost no 
restrictions in the transverse space [28]. Following the notation in eq. (17), an explicit 
form of the cut loop momentum is given by [28] , 

d—n+l d—n+l 

t^V^ a'nt, aW = -y^ (52) 

1=1 i=l 

where represents a specific linear combination of vectors {Ki\ and lies within physical 
space. The coordinates ai may take complex values and parametrize the transverse part 
of the loop momentum; only the square of the transverse part is fixed. We will denote 
the parameter space, a = {a;*}j=i,d_„+i with = —V"^, by in the following. It 
turns out that the remaining freedom in the transverse momentum directions allows us 
to uniquely specify all scalar and tensor coefficients in eqs.(26), (27) and (28). 

The numerator tensors contracted with the on-shell loop momenta (52) give linearly 
independent functions on the parameter space M^. This can be seen considering the 
definition of the basis of numerator tensors. For a non-vanishing value oi V"^ 
the parameter space contains a (d — n)-dimensional sphere, Sd-n C M^. The basis 
tensors (20) are symmetric and traceless in the transverse vectors n*'^. Thus, when we 
contracted a numerator tensor with the loop-momentum (52) we use ■ Kj = = -V 
and obtain harmonic polynomials in terms of «j. These polynomials are of course 
linearly independent being in one-to-one correspondence to spherical harmonics on the 
respective sphere, ^^-n- 

For example, for two-particle cuts in four dimensions, n = 2 and d — A, we may 
use spherical coordinates for the coordinate vector a in (52) and obtain the classical 
spherical harmonics YJ^, 

a = i\V\ X {cos^, sin ^ cos 0, sin sin 0} , (53) 
n'l^^ . . . nl^^^ . . T'" ~ Yi,^i9, 0), I = k ,m = -I, . . . ,1 , (54) 

as observed in [21] using a specific representation of the on-shell momenta. Definitions 
of curly-brackets and vectors in (54) are given in eqs. (17) and (20), respectively. The 
norm \V\ is given by \V\ = VV ■ V. Similarly, one obtains a representation in terms of 
spherical harmonics for generic dimensions and tensor-rank of the numerator tensors. 



Susy Theories and QCD: Numerical Approaches. 



27 



Quadruple cuts in four dimensions, n = d = 4, are a special case. For this case the 
transverse space is one-dimensional, «i = ±i|V^|, with the harmonics degenerating to 
two functions; the even and the odd functions on two points, ai = ±i\V\, representing 
a zero-dimensional sphere, Sq. The even and the odd functions, Y and Y'. on these two 
points can be defined to take the respective values, Y = {1, 1} and Y' = {1, —1}. The 
constant function Y corresponds to the scalar integral and Y' to the hnear tensor in 
eq. (26). An analogous behaviour appears as well away from four dimensions, d ^ 4, 
whenever we have d — n. 

A refinement of the parametrization in terms of spherical coordinates (54) is 
necessary when the square of vanishes; \V\ = 0. This case arises for cuts of integrals 
with massless internal and at least one massless external leg. The parameter space then 
takes the form of a cone, 

d-n+l 

^ q;W = 0. (55) 

i=l 

— * — * 

A generic parametrization of (55) is given hy a — {t, i tp} with t and (3 complex valued 
and = 1. Also in this case the numerator tensors have a unique functional form 
in terms of t and spherical harmonics originating from angular dependence oi /3. In 
particular, they are linearly independent functions on the parameter space M^. 

A special case worth mentioning are triple cuts with massless external legs; d = 4 
and n = 3. For this situation we have only a two dimensional transverse space, with 
i — 1,2, and the angular dependence through /3 degenerates further to /3i — ±1. The 
solution space to eq. (55) then consists of two branches, a — {t, ±it}. A simple example 
of this situation is given in section 4.4. Both branches of phase-space have to be taken 
into account to solve the unitarity relations. 

In summary, we see that the tensor numerators in eq. (11) are in one-to-one 
correspondence with a basic set of functions (e.g. spherical harmonics) that appear 
in the unitarity cut. In particular, all tensor coefficients may be identified using the 
linear independence of these functions. To this end, we may evaluate the unitarity 
relations (51) on a predefined set of points on the parameter spaces Mg and use the 
known functional form of each of the numerator tensors to pick out their coefficients. 
The degrees of freedom left after the on-shell conditions are imposed, are thus sufficient 
to uniquely obtain the loop integrand through unitarity cuts. 

4-3.1. Power Counting. A further interesting observation is, that the total angular 
momentum quantum number that appears, i.e. / of Yj, is constrained by the power- 
counting of the theory. For terms with P powers of loop momentum we have, 

Yi: I < P for maximal power of {l^)^ . (56) 

Conversely, the maximal angular momentum reflects the UV behaviour of the respective 
amplitude. For the special case oi d = n this reasoning does not apply. For simplicity, 
we focused only on the situation with |V^| 7^ and we only assert that a similar counting 
argument can be extended to the case \V\ =0. 
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(a) (b) 

Figure 6. The triangle integral function Z3(fc2, fcs, i^45i) (left) and the associated 
triple cut (right). The virtual states in the loop are fixed to be a complex scalar. 
While the scalar integral coefiicient turns out to vanish, some tensor coefficients are 
non-zero. The tensorial coefficients are important input for the computation of lower- 
point integral coefficients in the numerical unitarity approach. 



The total number of numerator tensors can then be counted using representation 
theory of the orthogonal group, SO{d—n+l), and are given by the number of harmonics 
up to a given total angular momentum. 

4.4- A. Triple- Cut Example 

We consider the amplitude ^45(1", 2~, 3^, 4+, S"*") with a scalar circulating in the loop. 
Amplitudes of this kind were first computed in refs. [22] using the unitarity method. 
Here we are only interested in the computation of a single triple-cut, as shown in fig. 6. 
It turns out that some tensor coefficients (26) are non- vanishing while the scalar triangle 
coefficient C234 vanishes. This simple example allows to illustrate how to apply some of 
the methods that were introduced earlier. 

In later parts of this review we will discuss further important aspects of such 
unitarity computations. In section 4.5 we will discuss, based on the present example, 
how to compute tensor coefficients numerically using discrete Fourier transformation, 
following the notation in [21]. Building on this, in section 4.6, we will explain how 
numerical precision can be monitored during computation of such triple-cuts. 

We start with the computation of the triple cut fig. 6. A momentum parametrisation 
and reference momenta, similar to eq. (17), are given by, 

n+ = A2A3 , = A3A2 , ■ Ki = , 

q = {X2 + tX,)~X2, r, = tXsX2, ^2 = A3 (-A3 + tA2) . (57) 

This is a special case of solutions of the on-shell conditions for integrals with massless 
corners as explained below eq. (54). The on-shell momentum space takes the form of 
a light cone and consists of two parts, with the second branch, related to (57) by 
parity reflection A* o A*. The loop momentum is then, 

= A2(A2 + t'A3) ^^ = t'A2A3, il = {-Xs + t%)Xs. (58) 
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Notice that a more standard way to write the above parametrization would use 



1 _ + n-)/^/s^, = {n+ - and if = -t^^{in^ ± n^). The 



n — i[n 

use of complex momenta is crucial in order to evaluate such unitarity cuts. 
The triple-cut, fig. 6, is given by the product of three tree amplitudes, 

c234(^o = ((-^2)^ 2-, 3+, ^^^r 4+, 5+, i-,m 

Inserting explicit tree amplitudes we then obtain, 
"'^'^'^^ [34] [^4 (-4)] [(-4) 3] 



{{-U) 4) (4 5) (5 1) (1 £2) {£2 i-U)) {{-£2) 2) (2 4) (4 {-£2)) 
and after using the above momentum parametrisation, eq. (57), 

, _ _ t^(13)^[3 2]((21)+t(31)) 

'''^ ^^ (3 4) (4 5) (51) 

C234(A')=0, (61) 

where the triple-cut vanishes on the second branch £^ due to the vanishing of the three- 
point vertices for these momenta. 

We can make several observations: (a) There is no term constant in t consistent 
with the absence of scalar triangles [22] . (b) The maximal power of t is three, consistent 
with the three powers of loop momentum expected from power-counting of a triangle 
diagram, (c) The expression (61) has no poles at finite non-zero values of t consistent 
with the absence of scalar boxes for this particular combination of external helicities 
and scalar states in the loop [22]. The equations (51) thus simplify with box-subtraction 
terms being zero and, thus, absent. 

The parametrization of the numerator tensors (3.3.3) is given by, 

C234(£) = ^ (n+ ■ £f + (n+ ■ £f + (n+ • £) + c° + 

c-^{n- ■ I) + c-2(n- • £f + c'^in- ■ £f . (62) 

The next step is to relate this parametrization of the numerator tensors (62) to the 
expressions for the triple cut eq. (61). For on-shell momenta we find for (62), 

C234(r) = (-S23t)' + ("523^)' + c' {s^^t) + C° , 

C234(^') = C° + {-S2,t'f + {-S2^t'f + C"^ (-^23^')' , (63) 

where any of the on-shell momenta £1'^ i = 2, 3 or 4, can be used. 

Comparing the polynomials in t and t' given by cut (61) and ansatz (63) we find 
all tensor coefficients, 

cO = = = c-S c^^2 = -S23 Af- ( ' . (64) 



(21) 

In an analytic approach, it is of course no problem to compare the t-dependent 
expressions. What matters here, is that the tensor basis gives linearly independent 
functions on the solution space of £"-''^, so that all tensor coefficients can be identified. 
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Figure 7. The discrete Fourier transformation to obtain integral coefficients. Tlie 
triple cut contribution depends only on a small number of functions. Instead of a 
contour integral, the triple-cut may thus be evaluated on a discrete set of points. All 
scalar and tensor coefficients are obtained exactly in a numerical procedure. 



This will turn out sufficient for the related numerical approach, as we will discuss 
in section 4.5. 

Thus, effectively by replacing t — )■ — (n"*" ■ l2)/s2'i in (61) we obtain the off-shell 
loop-momentum dependence, 

r (P\ . .tree ^(!^^:1MV ^ ^ ll^lllM^ ^fi.^ 

c,,m = -S2,A, ^21) ) (21) ) ' ^^^^ 

For the special case of this example, many coefficients including the scalar coefficient 
c° vanish. Because of this the associated triangle integral does not contribute to the given 
amplitude. The computation of the tensor coefficients is important, nevertheless, as they 
are used when computing the two-particle cuts, for which they serve as subtraction terms 
as apparent from eq. (51). 

In the above analysis we relied on the fact that we had analytic expressions available. 
In the following section we will see that a purely numerical approach can be set up to 
follow very similar computational steps to obtain scalar and tensor coefficients. 



4-.5. Discrete Fourier Transform 

We continue the discussion of the above example in section 4.4. A direct way to extract 
the tensor coefficients from the cut expression C234(^") in eq. (61) is through contour 
integrals, 

c'= = ^/yf^V'c234(t), (66) 
2-^1 J t \-s23J 

which take the form of a Fourier transformation in an angle for an integration contour 
t = toe'^. 

From a power-counting argument we know a priory that only a finite number of 
monomials in t may appear. In particular, we know that t will at most appear as a third 
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power t^. We may thus use a discrete version of the above contour integral fig. 7, 



giving an exact numerical way to obtain the required tensor coefficients. Typically the 
expected tensor-rank determines the number of points that need to be sampled in the 
above sum. For the example in section 4.4 four points, m — A would be sufficient. 

For the computation of all tensor coefficients in equations (26), (27) and (28) as well 
as in equations (88), (89), (90) and (91) analogous discrete sums can be constructed [21]. 
The central observation is that the numerator tensors are related to a finite number of 
functions. For phase-space integrals the orthogonality of these functions can be used to 
extract a particular tensor coefficient from a unitarity cut. For numerical approaches a 
related discrete versions of such integrals can be setup. 

4-6. Numerical Stability 

A clear understanding of the relation between numerator tensors and a function basis 
in phase-space can be further exploited; we can check the numerical precision during 
the computation of loop amplitudes. 

Loss of precision may occur on a given phase-space point, for example, due to 
an unfortunate choice of reference momenta (17) or a vanishing Gram determinant 
(35). A particularly universal way to deal with loss of precision is the use of higher 
precision arithmetics [79]. As demonstrated in ref. [21], for an efficient approach one 
tries to identify the computational steps that lead to large round-off errors and then 
to re-evaluates only the problematic contributions to the amplitude (and only those 
terms) using higher-precision arithmetic. Use of higher precision arithmetics is more 
time consuming. Such an approach requires that results be sufficiently stable in the 
first place, so that the use of higher precision is infrequent enough to incur only a 
modest increase in the overall evaluation time; this is indeed the case. 

The simplest test of numerical stability [79, 28, 21, 68, 38] is checking whether the 
known infrared singularity of a given matrix element has been reproduced correctly. 
Typically, a combination of various checks is required. A refined test [21] can be 
setup to check the accuracy of the vanishing of certain higher-rank tensor coefficients. 
From power-counting arguments we know on general grounds which high-rank tensor 
coefficients have to vanish. All tensors with rank greater than m must vanish, for the m- 
point integrals with m = 2, 3, 4. The central advantage of this check is that small parts 
of the computation may be singled out for re-computation, reducing the computation 
time accordingly. 

For the above example of a triple cut computation, section 4.4, a check of the 
vanishing tensor coefficients may be implemented as an extension of the discrete Fourier 




with m > 4 , 



(67) 
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Figure 8. The box integral function I{Ki2) (left) and the associated quadruple cut 
(right). The helicity sum over internal states is suppressed in the graphic. 

sum. If a tensor of rank-four were to exist it could be computed numerically by, 

1 ™ /j- 2TTij/m \ 

^fc=4 =±yl C234(t = toe'"'^/'") , with m > 5. (68) 

Clearly should always turn out to be zero, c'* = 0. For finite precision the deviation of 
from zero can be used to monitor precision of the intermediate computational steps. 
This test requires to sample more values of the complex parameter t; m = 5 instead of 
m = 4. However, observing an instability at the level of an integral coefficient as opposed 
to a failing check of an infrared-pole has the big advantage, that one can recompute very 
targeted a small part (the specific triangle coefficient) of the full amplitude. 

The parameters of such a rescue system have to be tuned in order to optimize 
efficiency. This includes the question which integral coefficients to check and what 
deviations from zero to expect for a given tensor coefficient. Some further details of this 
approach can be found, for example, in [21]. 



4-.7. A Box Example. 

Here we will discuss the pure gluon amplitude y4|(l^, 2^, 3+, 4+, 5+). The analysis below 
will be purely analytic, however, in the way we set it up it is straightforward perform 
all computations numerically. Since this amplitude has five external massless legs, the 
box functions can have one external massive leg, and three massless ones. One of these 
box functions with massive legs K12 = ki + /c2, is given by. 



d^-'H 1 



(27r)4-2^ £2(£ ^2)2(£ _ _ ^2 - + ^5)2) 



and displayed in fig. 8. The other boxes have the massive leg being -^34? -^45 o^' 

^51- 

We focus here on the analysis of the coefficient of X{Ki2). The first step is to 
parametrize the integrand (ii345(£) in eq. (11) that extends the box function with a 
numerator tensor. The transverse momentum space (17) is one-dimensional and is 
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spanned by a single vector n^^^, 

■Ki2 = n^ ■h = n^ -h^n^ -h^O, (69) 

with the exphcit solution for n^'^ given by, 

n'-£ = ([3\hhe\3) - (3| W|3])/s34S45 , (70) 

where the loop momentum i is identified with £i as defined in fig. 8. The only symmetric 
traceless tensor (20) in the transverse space is given by the vector n^'* itself; — n^^ 
as in section 3.3.3. The numerator tensors of the box integrals are then parametrized 
as, 

resembling the generic form in eq. (26). 

The coefficients ^5^345 and ^^345 have to be determined by comparing the quadruple 
cut equation (51) to the ansatz (71) at the associated on-shell kinematics. In four 
dimensions the on-shell conditions, 

^2 ^ = (£ - ^5)' = (^ + h^f = (^ + K^A? , (72) 

have two solutions. A general form for the on-shcll momenta can be found in ref. [21]. 
Explicitly, one set of on-shell momenta is given by, 

solution (a): 4 = ^^^^A^A^ , 4 = -^^|^A,(^,,A3) , 

with the second set of on-slicll momenta, solution (b), being related by parity 
conjugation; exchanging Aj •H- Aj. We used here notation as discussed in eq. (32). 

The quadruple cut is given by a product of four tree amplitudes summed over 
internal helicity states, 

diM,{Q = J]Af^((-£2)-"^l-,2-,£3'^)Ar((-4)-'^3+,£^) X 

hi 

X Af^%{-h)-''\A+,£'l')Af'^{{-£i)-''\5+,£'^') . (74) 

Only the momenta in eq. (73) give a non-vanishing contribution to the quadruple-cut 
(74); for the internal helicities {hi, /12, h^, h^} = {—, — , -|-, -|-}. In total we have, 

solution (a): J,34.«.) = (24) (4 (-«>{(-<.) 1) [<4 (-Im-^) 3] ' 



(4£i) ((-£4)4) [5 £2] [£2 (-4)] 
= ^S34S45^r(l",2-,3+,4+,5+), 
solution (b): (ii345(£j) = . (75) 
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In order to determine the integral coefficients 0^1345, d\^^^ we have to evaluate ansatz 
for the integral on the unitarity cut. For the particular momenta no other terms in the 
ansatz contribute leading terms in the factorization limit. We thus obtain, 

rf"l345(C'') = ^45 + ^1345 ' ^f) = ^45 ± ^^345 , (76) 

where any on-shell momentum £j may be used in the numerator tensor, due to the 
properties of the vector v}'^ os specified in eq. (70). 

Using the values of the quadruple cut in eq. (75) we thus can solve for the unknowns 

^1345 ^1345' 

rf?345 = ^1345 = ^^34^45 ^f^ 2', 3+, 4+, 5+) . (77) 

That the integral coefficients ^^345 and d\^^^ turn out to take identical values is a low 
point accident related to the presence of three-point amplitudes in the unitarity cut. For 
the value of the one-loop amplitude the tensorial term is of no immediate importance; 
it drops out after integration. However, for the computation of the remaining triangle 
and bubble coefficients this is an important ingredient. 

The normahzation factor 1/2 in scalar integral coefficient (77) is usually attributed 
to an 'averaging over' quadruple cut solutions [25]. In the presented computation this 
factor appears, somewhat differently, from inverting the set of linear equations (76). 

4-8. Rational Terms From D-dimensional Unitarity Cuts 

The rational terms (7) are free of branch cut singularities and cannot be detected 
by unitarity methods in four dimensions. Away from four dimensions, D = 4 — 2e, 
rational terms carry factors of (— s)~^ = (A)^^(l ~ eln(— s//t^) -|- 0{e^)), in order 
to compensate for the dimensionality of the coupling constant. The so introduced 
logarithms, in turn, make it possible to detect rational terms via D-dimensional unitarity 
methods [23] (see also the early work [129]). This version of unitarity, in which tree 
amplitudes are evaluated in D dimensions, has been used in various analytic [72, 130] 
and numerical [29, 131, 124, 132, 79, 68, 133, 39, 70] studies. For a detailed discussion 
of analytic D-dimensional approaches we refer to [43] . Here we will follow the approach 
discussed in [29] including elements of [124]. 

4.. 8.1. Dimension Dependence. The prescription of the D-dimensional unitarity 
method is to compute the D-dimensional loop-amplitudes and in the end take the limit, 
D = 4 — 2e with e — > 0. The combined dependence of integrals and their coefficients on 
the regulator e yields the finite rational terms. 

In numerical approaches the dependence on the dimension parameter D is not 
as accessible as in an analytic approaches. Amplitudes may only be computed in 
a fixed dimension. The typical strategy, to deal with this issue, is to numerically 
compute in distinct discrete dimensions. With a clear understanding of the dimension 
dependence [29], the dimension parameter may be reinstated in the final steps of the 
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computation. The limit D — )■ 4 is then performed analytically. We will discuss the 
D-dependence in detail in the following. 

The origin of the D-dependence is twofold. While the external momenta and states 
are kept in four dimensions, polarization states and momenta in the loop are extended 
beyond four dimensions. The two sources for dependence on D are the sums over 
virtual polarization states in D dimensions and momentum invariants formed using 
D-dimensional loop momentum. 

We consider first the dependence of rational terms on D-dimensional polarization 
states. A priori one would need to sum over the particle spectra of virtual gluons and 
fermions extended beyond four dimension. This can be done in practice [29], however, 
we will describe a shortcut, that avoids considering polarizations states beyond four 
dimensions altogether. Applying the arguments of section 3.5 we can relate the rational 
terms of QCD amplitudes with internal gluons and fermions to amplitudes with virtual 
scalar and fermionic states. The extension of scalar states to D-dimensions is straight 
forward as, in contrast to gluonic states, no additional polarizations have to be taken 
into account. In fact, it turns out that the D dependence originating in the change 
of the number of polarization states with dimensionality my thus be avoided. Such an 
approach is computationally more efficient. The computation time grows faster than 
linearly with the number of particle states circulating in the loop. The computation 
of amplitudes of a complex scalar as opposed to (massless) vector particles leads to 
efficiency gains in D > 4 dimensions. 

For the computation of rational terms it is, thus, sufficient to consider simplified, 
yet generic amplitudes with: 

(i) closed scalar loops, 

(ii) mixed scalar and fermion loops. 

Even though QCD has no scalars, introducing them is a useful trick. For QCD like 
theories these assumptions are no restriction for the computation of rational terms; 
contributions from either only virtual gluons or fermions are related to virtual scalars 
(44). Similarly, rational terms from mixed gluon and fermion loops can be mapped to 
the ones with the gluon replaced by a scalar, see e.g. [78, 24]. Loops with mixed fermion 
and scalar states do only give rise to state sums over fermions which exit to external 
sources. Since for these amplitudes no closed loop of Dirac gamma matrices can be 
formed, no explicit dependence on the number of fermion states appears. We thus avoid 
the dependence on the dimensionality through the number of spin states. 

The second source of the dimension dependence of loop amplitudes is the 
dependence on the dimensional loop momentum. This dependence is simplified due 
to the fact that we keep external polarizations and momenta in strictly four dimensions. 
Splitting the loop momentum into a four- and {D — 4)-dimensional part, — -\- 
the dependence on the {D — 4)-dimensional part is limited to the form. 




(78) 
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That is, no other vectors with non- vanishing components in — 4) are available to form 
scalar products and we ignore linear /x-terms built from e-tensors. Rotation symmetry 
in the {D — 4)-dimensional part of momentum space is preserved in this way. 



4-8.2. The Loop Integrand. The integrand basis has to include pentagon terms, when 
working in D-dimensions 



(^) 

W - 2^ r., r.. r,. r,. D. + 2^ 



h<...<h -^w A2 A3A4A5 ii<j2<i3<i4 Di^Di^Di^Di^ 

+ 2^ D,,D,,D,, + ^ D-D- ' 

with propagators and numerators depending on the £)-dimensional loop momentum 
_ £^l £n_ suppressed the tadpole terms ai^(^) which are not needed for 

computation of rational terms in amplitudes with only massless states. As compared to 
the earlier expression (11) the form of the numerator terms has to be adapted to the 
D-dimensional will be discussed in section 4.8.4. The generalization to massive 

states has been presented in [132]. 

The absence of terms with more than five propagators in eq. (79) is due 
to the restriction to strictly four dimensional external momenta and polarization 
vectors [115, 116]. (See also [117] for a recent discussion including generalizations to 
two- loop integrals.) Integrands with more than five propagators can be reduced to at 
least pentagon terms. At most five independent momentum vectors can be formed in 
D dimensions, since the (D — 4)-dimensional momentum is conserved in the loop. 
Because of this, Gram determinants which depend on the loop momentum and five or 
more independent momenta have to vanish. The resulting identities between momentum 
invariants can be used to reduce higher n-point propagator structures recursively to at 
least pentagon integrals. 



4.8.3. D -dimensional Unitarity Relations. As for the cut-containing parts, the 
computation of D-dimensional loop amplitudes amounts to determining the coefficients 
of a basis of scalar and tensor integrands and subsequently performing the loop 
integrations. A generalization of the unitarity relations (51) can be used to completely 
determine the loop integrand (79). 

The restriction to purely four-dimensional external momenta has important 
implication for the unitarity cuts of the internal propagators. In fact, the {D — 4)- 
dimensional component of the loop momentum is conserved and enters all propagators 
in the form of a mass term. 



For unitarity cuts, this observation leads to an important restriction. In addition to the 
four momentum components of we thus obtain only one additional degree of freedom. 
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i.e. /i^. We can then put at most five propagators on shell, allowing us to consider at 
most penta cuts. Of course this meshes well with the absence of higher point functions 
in eq. (79) in the first place. 

The unitarity relations are generalized to D dimensions, by including an additional 
generalized penta cut, 

eni.isi.iM = E <r(^)<^W<^(^)<^W<r(^) ' (81) 

D-dim states 

which is defined for appropriate on-shell loop momentum, £. The remaining unitarity 
cuts are implemented similar to the four-dimensional case (sec cq. (51)), however, 
including additional subtraction terms from the pentagon-coefficients GiiiihiiirX^) ■ 

We may work in any discrete dimension bigger than four, which accounts for the 
full loop-momentum dependence in D dimensions. The internal state sums have to be 
extended to the D dimensional setup. For the case of a scalar state circulating in the 
loop, the D-dimensional extension of the scalar tree amplitudes is to be considered. 
Given the simple dependence (78) on the {D — 4)-dimensional loop momentum it 
is sufficient to consider a scalar in five dimensions (see also [29]) to obtain the full 
dependence on yu^. Equivalently, one can consider four-dimensional virtual scalar states 
with a (dynamical) mass [23, 124]. 

Typically, the on-shcll conditions do not constrain the loop momenta completely, 
but lead to a variety of solutions. Enforcing the unitarity relations (81) on the variety of 
on-shell solutions gives an infinite set of equations. We will discuss the parametrization 
of the D-dimensional numerator tensors below in section 4.8.4 and section 4.8.5 which 
allows to determine all needed integral coefficients. 



4.8.4- Numerator Tensors. In order to parametrize the numerator tensors of a given 
propagator structure in eq. (79) it is convenient to introduce an adapted vector basis 
in momentum space. The D-dimensional momentum space is split into three subspaces 
spanned by the vectors Xj, and [29], 

Ki for i G { 1 , ■ • • , n — 1 } , 

n' for « G {l,---,4- (n- 1)} m' for ie{5,---,D} (82) 
{n\ n^) = 6'^ , {m\ m') = S'^ , {n\ m?) = , {n\ Kj) = , 

4-(n+l) D 

9T^ E g^n^'-^^T.^^V- (83) 

i=l i=5 

In a first step, the £)-dimensional space is split into two subspaces: the 4-dimensional 
subspacc and its {D — 4)-dimensional complement, m*. In a second step, the 4- 
dimensional momentum space is further decomposed into a physical space, spanned by 
the independent external momenta of the integral Ki and their transverse space within 
four dimensions, {n*}. For convenience the bases {n*} and {m*} are orthonormal. The 
metrics g'^ and g'^_4 are projections of gj^i, to the respective subspaces. 
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We can follow similar steps as in section 3.3.2 to conclude that the generic numerator 
tensors are traceless symmetric tensors in the whole of the transverse space, 

where the curly-brackets indicate symmetrization and subtraction of all traces. More 
explicitly, the traces should be subtracted in transverse space, using the metrics g'^ and 
fl'D-4' 9^" ■ rotation invariance in (D — 4) dimensions, the dependence 

on the m* is further restricted to appear solely in terms of g^_4^ = '^f=5 m^^m^" . 

A subtlety appears here for numerical applications; the traceless- condition 
introduced a dependence on the dimensionality of space time, as illustrated by the 
following example, 

D DA 
n^,u = ^M^^^} = {i9D-4)^u - ^_ . ^ {9±)^cu) , (85) 

i=5 ^ ' 

yielding a D-dependent coefficient of g^^ and thus an explicit D-dependence. For 
numerical computations we would like to avoid specifying D and keep it as a parameter. 
Without loss of generality, the dependence on the space-time dimensionality can be 
tracked more conveniently by allowing terms with non-vanishing traces. This can 
be achieved at the minor cost of introducing terms that do not integrate to zero, 
section 4.8.6. In particular, as given in [29] we use instead a representation of the 
numerator tensors in terms of tensors that are trace-free only in the physical transverse 
space, 

'Vl - Mfc — "{fji ■ ■ ■ ■ ■ ■ 9fik-itik ' 'VlM2-/ifcy± — U ) 

The above rank-two tensor (85) would then be given by, 

n^^rr = {gD-4)^^J''e'' = , (87) 
without explicit dependence on the parameter D. 

4-8.5. Tensor Basis. An explicit form of the numerator tensors was given in [29]. The 
pentagon numerator tensors are, 

where the absence of any loop-momentum dependence has to be noted. For this case the 
transverse space coincides with the {D — 4) dimensional space and no traceless tensor 
can be formed. Therefore terms of the form /x^ can be converted into inverse propagators 
and scalar terms and arc represented by lower-point integrals and the scalar pentagon. || 
The remaining numerator tensors are, for the box, 

II Ignoring e-tensors ^-dependent terms of the pentagon can be converted to propagators and scalars; 
— /i^ = igD-^^tivi^t" — — Ylii=\ ii^i ■ ^)(^* ■ ^)- Terms of the form {Ki ■ I) can be further expressed 
in terms of inverse propagators as in eq. (22). 
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for the triangles, 

Cni^is (i) = Qu.i3 (^) + {cUish + cl,^i,t2 + cl,^i,) , (90) 
and finally for bubble coefficients they are, 

~k,,,ii) = KM+f,%l,^. (91) 

Here the vectors we introduced = (n* ■ i). The are defined for each propagator 
structure of (79) individually as defined in eq. (83). The tadpole contributions may 
be found in the original literature. The /^-independent tensors take the same form 
as given earlier, being traceless symmetric tensors in the vectors rf. These tensorial 
structures parametrize the most generic tensor integrals, subject only to power-counting 
requirements of QCD like theories as well as the constraints from rotational symmetry 
manifest in dimensional regularization. 

4-8.6. Integration. The final step to extract the D-dimensional amplitude is to evaluate 
the integrals. We are interested here only in the computation of the rational terms. Only 
the new integrand structures of eqs. (88), (89), (90) and (91), proportional to powers of 
jj, may contribute to the rational remainder in eq. (8). 

As above terms dependent on vectors integrate to zero due to their angular 
dependence and can be dropped. In order to obtain the rational contributions we are 
left with terms proportional only to powers of fj?. Without speUing out the details, the 
non- vanishing hmits are given by [134, 23, 124], 

, d^i 1 /• d^l fj,^ 1 

lim / — ^-7- = , lim — 



J in'^/^ Di,Di,Di,Di, 6' J in^/^ Di,Di,Di, 2' 

hm / — FTTTT— ^- = Kf^, (92) 

The limits of these integrals combined with the associated integral coefficients add 
up to the rational term of the one-loop amplitudes. The rational term is then given 
by [130, 131, 29], 

ii<i2<i3<i4 ii<i2<i3 ii<i2 

We note that the D-dimensional approach is very general and can be applied to 
higher loop computations. For evaluation of scheme shift and order e contributions we 
refer to [29]. 



5. On-shell Recursion 



On-shell recursions [47] rely on on-shell scattering amplitudes with a fixed number of 
partons in order to obtain the ones with arbitrary multiplicity. The underlying structures 
used in this approach are universal factorization and analyticity properties. 

In fact, when intermediate states are nearly on-shell, amplitudes factorize into 
products of lower-point amplitudes. At tree-level, the naive attempt to invert 
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factorization equations and assemble the on-shell factorized amplitudes into their 
parent amplitude raises two central questions: Firstly, how to recover the full off- 
shell kinematics away from the factorization limit? And, secondly, how to combine 
various factorization channels while avoiding double-counting? On-shell recursion was 
introduced by Britto, Cachazo, Feng and Witten [47]. These naive obstructions are 
overcome by the effective use of complex kinematics and the use of Cauchy's residue 
theorem. 

While originally constructed for tree amplitudes on-shcll recursions can be extended 
to loop- level [19, 67, 20]. Compared to tree amplitudes, the obvious difficulty is that 
loop amplitudes contain branch cuts, which complicate the use of Cauchy's theorem. 
The resolution to these difficulties has been to work on subparts of the amplitudes, 
which are purely rational functions. The price to pay is that additional, un-physical 
singularities have to be dealt with, which are present in subparts but cancel once 
the full amplitude is assembled. Un-physical poles can typically be understood from 
complementary information. In the case of on-shell recursions for rational terms spurious 
singularities are characterized using prior knowledge of the logarithmic parts of the 
amplitude. In addition, while factorization of loop amplitudes is understood for real 
kinematics, presently there are no theorems on the factorization properties of loop 
amplitudes with complex momenta. Indeed, there is a class of poles, 'unreal poles' 
[19], whose contributions have to be taken into account and whose nature is not yet 
fully understood. For most applications unreal poles can be avoided using alternative 
factorization channels. (See, however, ref. [135] for recent progress with understanding 
the origin of unreal poles.) 

Most recent developments focus on recursion relations at the integrand level [136] 
and have already led to many results in maximally supersymmetric Yang-Mills theory. 
Here we are mainly interested in on-shell recursions for loop amplitudes, which can be 
applied after integration has been performed. 

For numerical applications at tree-level, on-shell recursions are efficient, allowing 
to repeatedly take advantage of the remarkable simplicity of the physical scattering 
amphtudes (e.g. Parke- Taylor amplitudes). At loop- level, on-shell recursions hold 
the potential of an efficient complement to unitarity approaches for two reasons: (a) 
Recursions can sidestep the computation of tensors coefficients (11) and deal entirely 
with the much smaller set of scalar integral coefficients, (b) On-shell recursions for 
rational terms rely on strictly four dimensional information and, thus, do not require 
more involved D-dimensional objects. Despite the potential of loop recursions, unitarity 
approaches are very universal and straightforward to implement for a large class of 
theories explaining their widespread use. The BLACKHAT-library [21], for example, 
makes use of both approaches and like this improves the efficiency of computations. 

We will first point out some central ideas of on-shell recursions at tree-level and 
will then turn to loop-level recursion. 
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5.1. Tree Recursion 

For a more complete discussion on tree level recursions we refer the reader to the 
chapter of the present review dedicated to this subject by Brandhuber, Spence and 
Travaglini [53] . Here we will merely introduce basic notation needed for the discussion 
of loop-level on-shell recursions. 

On-shell recursion relations systematically construct a scattering amplitude (or 
parts of it) from its poles and residues in momentum space. Typically, a parametrization 
of the complexified phase-space in terms of a complex parameter z is introduced. An 
amplitude An{z) = An{{ki{z)}) is then reconstructed as the unique function with 
consistent poles and residues in the 2;-plane. Analyticity is in many cases so restrictive for 
the form of the amplitudes, that this is in fact possible. Prior knowledge of singularities 
in phase-space as well as a simple physical interpretation is a necessary input in this 
procedure. 

A particularly useful parametrization of a plane in complexified phase-space is given 
in terms of a deformation of two momentum vectors [47] , 

An{z) = An{ki,...,kj{z),kj+i,...,ki{z),...,kn). (94) 

The momenta ki{z) are chosen to keep momenta on-shell and the overall momentum 
conservation intact. The explicit form of such a parametrization, denoted by is 
given by the linear transformation, 

[;', I) : Xj Xj - zXi , Xi^ Xi + zXj , (95) 

where 2; is a complex parameter. This shift of spinors leaves untouched A^, A;, and the 
spinors for all the other particles in the process. The corresponding momenta are, 

kj{z) — kj — z XjXi , ki{z) — ki + z XjXi , (96) 

and on-shell conditions, here k^{z) = = kf{z), as well as overall momentum 
conservation remain intact, due to kj{z) + ki{z) — kj + ki. 

Poles in the variable z appear through propagators when intermediate momenta go 
on-shell, 

-±- = i (97) 

where a denoted a range of momentum vectors that contain momentum ki but not kj. 
The location of the pole is given by z^ — —Kl^/ \l\Ka\j). Near the singularity of the 
propagator (97), the amplitude is given by its factorization properties, 

and uniquely defines the residue on the pole in terms of the on-shell lower-point 
amplitudes A^^(2;„) and A~\{za). 

By using Cauchy's residue theorem we can relate the amplitude at 2; = 0, A„(0), to 
its poles. For a vanishing contour integral that encloses the point z = 0, we find, 

O^i^^ ^ ^(0)=_ ^ Re.„,„^. (99) 
Jc 2m z ^-^ z 

poles a 
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The vanishing contour is associated to large complex momenta. In this region power- 
counting arguments can be used to show the vanishing of the function An{z — )■ oo) = 0. 
One may also include the residue at 2; = oo. However, this residue is not related to a 
factorization limit as above (98) and has to be obtained by other means (see e.g. [67]). 

When there is no large-2; contribution, the physical amplitude An{0) is obtained 
explicitly by, 

MO) = E E <L{Za)^A-%{z^) . (100) 

ah " 

The on-shell amplitudes with fewer legs, A'^ ^ and A~^j^, are evaluated in kinematics 
that have been shifted by eq. (95) with ki{za), such that their intermediate momentum 
is on-shell, Kaiza)"^ = where = —KfJ [l\Ka\j) ■ 

The n-point amplitude An is thus expressed in terms of sums over on-shell, but 
complex continued, amplitudes with fewer legs. These recursion relations can be 
extended to massive QCD, more general other theories and beyond four dimensional 
applications [137, 138]. Moreover, for certain hehcity configurations, they lead to new 
all-multiplicity expressions for these amplitudes [139]. 



5.2. Recursion for Loop Amplitudes 

On-shell recursion have been generalized to loop- level in refs. [19, 67, 20]. Following 
a similar strategy as outlined in section 5.1 for tree amplitudes, loop amplitudes are 
analyzed as functions of a complex parameter z. When applying a shift (95), we obtain. 



^l-loop/ 



Cniz) + Rniz) . (101) 



where Cn{z) denotes the logarithmic part of the amplitudes, and Rn{z) the rational 
remainders (8). In contrast to tree amplitudes, loop amplitudes (101) typically have 
branch cut singularities in the variable z originating from the logarithms within Cn{z). 

The way we deal with branch cut singularities is to focus on subparts of the 
amplitude which are free of logarithms: 

(i) the rational remainder: Rn{z) 

(ii) the integral coefficients: d^{z), c^{z) or b^{z) within Cn{z) 

(iii) D-dimensional coefficients: e^(^), d^{z), c^{z), b^{z) or a°(-z). 

These integral coefficients are introduced in sections 3.3.3 and 4.8.5, respectively. These 
terms can be shown to inherit factorization properties on physical poles from universal 
factorization relations at loop-level. For an intermediate propagator going on-shell there 
are typically three contributions for each internal state, 

lim /1 1- loop _ 4 tree ^ 4 1- loop , 4 1- loop ^ /i tree , 4 tree ^ 4 tree (^r)0\ 

In the first two terms, one of the factorized amplitudes is a one-loop amplitude and the 
other is a tree amplitude. The last term will appear here simply as a one-loop correction 
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Figure 9. Using Cauchy's ttieorem, rational expressions in loop amplitudes can be 
reconstructed from residues at poles in the complex plane. The poles are of two types: 
physical and spurious. All pole locations are known a priori. Residues at physical 
poles are obtained from universal factorization relations (102). Residues at spurious 
poles are obtained from the cut parts. 

Reprinted jig. 9 with permission from [21 j p. 16. Copyright (2008) by the American Physical Society. 



to the propagator. (For details about subtleties of the 'factorization function' in 
massless theories we refer to [140].) 

When considering subparts of the full amplitude, additional singularities, called 
spurious singularities, may appear. These singularities naturally cancel out in the full 
amplitude. A method to compute contributions from spurious singularities to rational 
terms R„ will be discussed in section 5.3.2. 



5.3. Recursions for the Rational Part. 

The terms Rn (8), which are purely rational in the spinor variables, cannot be computed 
using four- dimensional unitarity methods. On-shell recursion, however, allows us to 
construct these terms from purely four-dimensional data. Typically the information 
contained in the recursion relation has to be completed with properties of the cut part. 

The central input for this approach is the detailed understanding of Rn{z) as 
a function of 2; we need to know the location of poles, their residues and an 
integration contour, along which Rn{z) vanishes. Poles in the rational terms Rn{z) 
may be grouped into two classes as shown in fig. 9: physical and spurious. The 
physical poles are present in the full amplitude An{z), and correspond to genuine, 
physical factorization singularities. The spurious poles are not poles of An{z) and 
cancel between the cut parts C„(z) and rational remainders Rn{z). They originate 
from the presence of tensor integrals in the underlying field-theory representation of 
the amplitude, and appear as Gram determinant denominators (35). In unitarity 
approaches, inverse Gram determinants enter integral coefficients through on-shell loop- 
momentum parametrization. The exponent of the inverse Gram determinant in integral 
coefficients can be bounded by power-counting arguments. These denominators give 
rise to spurious singularities in individual terms. 

An example for the appearance of Gram determinants in unitarity approach can 
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be seen in the quadruple-cut computation in section 4.7. The Gram determinant (see 
also [21]), 

A(Xi2, h, h, h) = 2s45 (4|i^i2|5] [4|Ki2|5) , (103) 

appears in the on-shell loop momenta (73), through the factor 1/ [4|i^ri2|5); the remaining 
invariants cancel. The second factor, 1/ [51X1214), appears in the parity conjugate on- 
shell momenta. This 'chiral' dependence on inverse Gram determinants is typical for 
unitarity cuts with massless corners. In this particular example, the dependence on 
the Gram determinant (103) cancels in the final form of the integral coefficients in 
eq. (77). For the hclicity configurations discussed the integral coefficients arc the same 
as in maximally supcrsymmctric Yang-Mills and thus are effectively supersymmctric. 
Supersymmetric cancellations then reduce the naive power-count of loop momenta from 
four to zero, such that no dependence on Gram determinant remains in the integral 
coefficient. Explicit dependence of box-integrals on Gram determinant singularities can 
be found for example in analytic expressions in ref. [141]. 

Separating the different contributions of the shifted amplitude (101), we may write, 

R^iz) = R^{z) + Rl{z) + R':^^^ \z) , (104) 

where Rn{z) contains all contributions from physical poles, Rn{z) the contributions 
from spurious poles, and R^^^^'^ ^{z) the possible contributions from large deformation 
parameter z, if Rn{z) does not vanish there. More explicitly, the rational terms can be 
expressed in terms of a partial fraction decomposition in 

a " (7=0 

«.n.)=i:(^-r#^). (-) 

where the coefficients Bp, C/j, D^j are functions of the external momenta. The poles 
in z in cq. (105) are shown in fig. 9. The physical poles labeled by a arc gcncrically 
single poles. ^ In general, in a rcnormalizable gauge theory, the spurious poles, labeled 
by /3, may be either single or double poles. If Rn{z) vanishes for large z, the D„ are all 
zero. If not, then Dq gives a contribution to the physical rational terms, i?n(0). 

5.3.1. Physical Residues. The contributions of the physical poles may be obtained 
efficiently using the on-shell recursive terms represented by the diagrams in fig. 10. The 
'vertices' labeled by 'tree' denote tree-level on-shell amplitudes A^^, while the loop 
vertices Hoop' are the rational remainders of on-shell (lower-point) one-loop amplitudes 
-Rm, m < n, as defined in eq. (7). The contribution in fig. 10(c) involves the rational 
part of the additional factorization function [140]. It only appears in multi-particle 
channels, and only if the tree amplitude contains a pole in that channel. Each diagram is 

*{ Some shift choices may lead to double poles [142]; we can generally avoid such shifts [67]. A different 
approach based on unitarity cuts was recently suggested [135]. 
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(a) (b) (c) 

Figure 10. Diagrammatic contributions to on-shell recursion at one-loop for a [j, /) 
shift. The labels 'Zoop' and 'free' refer respectively to (lower-point) tree amplitudes 
^tree rational remainders of one-loop amplitudes R. The central blob in (c) is the 
rational part of a one-loop factorization function F [140]. 



associated with a physical pole (97) in the z plane, as in fig. 9, and is computed similarly 
to the recursive diagrams at tree-level (100). Further details for the computation of the 
recursive diagrams in fig. 10 has been described in refs. [19, 143, 111]. 

5.3.2. Spurious Residues. One approach to compute spurious residues is to use a 
completion of the integral functions, so called 'cut-completion' [19, 67]. To this end, the 
integral basis is adjusted to subtract off spurious poles within their integral coefficients. 
This in turn moves all spurious poles from the rational part i?„ in (101) to the cut 
part C„, such that the redefined rational remainder has poles only at the position of 
physical factorization singularities. The coefficients Bp and vanish in this approach. 
Attention has to be paid to the factorization equation which will mix cut part and 
rational part non-trivially. This approach has led to the computation of the rational 
terms for a variety of one-loop MHV amplitudes with an arbitrary number of external 
legs [19, 67], as well as for six-point amplitudes. 

For the purposes of a numerical program, however, it is simpler to extract the 
spurious residues from the known cut parts [21]. These residues, being absent in the 
full loop amplitude, are guaranteed to be the negatives of the spurious-pole residues in 
the rational remainder. That is, the spurious contributions are, 

hm Cn{z) = - f — ^ + + 0{z) , (106) 



z^ZB \{z — Zfs)'^ Z — Zp 

where Cn{z) is the shifted cut part appearing in eq. (7). Terms containing logarithms in 
the kinematic invariants cancel, but may as well be ignored for simplicity. The spurious 
poles at Zjs correspond to the vanishing of shifted Gram determinants, Am(2;^) = for 
m = 3,4, associated with triangle and box integrals. 

In order to compute the rational parts Rn we need to extract all residues Bp and 
Cp from the cut pieces Cn{z). For this we evaluate the integral coefficients di{z),Cj{z) 
and hk{z) numerically for complex, shifted momenta in the vicinity of the spurious pole. 
We also need to evaluate the loop integrals. This is done after using as input analytic 
series expansions of the integrals around vanishing Gram determinants, 

Xi{z) — ^ + logarithms 

Xf(^) = {z-zpf[p, + {z-zp)p2 + ..). (107) 
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Figure 11. We obtain the residue at the spurious pole located &t z — zp in the 
complex z plane by a discrete Fourier sum, evaluating C^{z) on the (blue) squares on 
the circle of radius centered on zp. The locations of other poles are represented by 
(red) dots. We ensure that Sp is sufficiently small so that other poles give a negligible 
contribution to the residue. 



The coefficients pi are rational functions in the kinematic variables which take a universal 
form. Examples of this procedure may be found in [21] with further details for generating 
these expansions using a dimension- shifting formula [144] may be found in [145]. It 
is important to have a precomputed analytic expansion available, as the logarithmic 
terms have to be dropped by hand; these terms do not cancel spurious poles in 
rational terms, but rather cancel between different integral functions. Thus we may 
avoid computing any logarithms or polylogarithms at complex momentum values. The 
expression obtained by replacing Cn{z) according to these rules, in the vicinity of z^, 
will be denoted by C^{z). 

5.3.3. Discrete Fourier Sum for Spurious Residues. The spurious-pole residues can be 
extracted from the cut parts by using a discrete Fourier sum. We evaluate C^{z) at m 
points equally spaced around a circle of radius 5^ in the z plane, centered on the pole 
location Zj^, as depicted in fig. 11; i.e., z = zp + 5pe^'^^^/^, for j = 1, 2, . . . , m. 
We can extract the coefficients Bp and Cp in eq. (105) via. 



Both m, the number of evaluations, and 6p are adjusted to optimize efficiency and 
numerical precisions. In general, an increase in m increases the precision, but at the 
cost of computation time. We choose Sj^ to be much smaller than the distance to nearby 
poles, but not so small as to lose numerical precision. 



Reprinted fig. 11 with permission from [21] p. 22. Copyright (2008) by the American Physical Society. 
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5.4- Recursion Relations for Integral Coefficients 

Here we address the question of whether we can also apply recursion relations to the cut 
containing pieces Cn{z). In fact it has been shown in [20] that it is possible to construct 
recursions, not for the full amplitudes, but in certain cases for the rational coefficients 
of the integrals. Such recursions can have important implication for the computation of 
rational terms, as we will discuss below in section 5.5. 

With unitarity methods available, one of the central motivations for coefficient 
recursion are efficiency gains and insights for the recursive computation of rational 
terms. Rational recursion relies on the data from cut containing terms; the residues of 
Gram determinant poles. Gram determinant singularities are tied to their associated 
integral coefficients. A direct recursive approach to compute integral coefficients directly 
gives access to spurious residues of Gram determinant poles. 

A firm grasp of the analytic properties of the integral coefficients is required. 
These properties differ from those of full amplitudes; in general these coefficients, 
being subparts of the full amplitude, contain physical poles as well as spurious poles 
of their own. Very much like for the case of rational terms, spurious poles originate 
from reduction of higher point tensor integrals and are associated to Gram determinant 
denominators. 

The spurious poles are harder to deal with in a purely recursive way. An intricate 
web of cancellations between distinct integral coefficients guided by the integral functions 
leads to the cancellation of spurious poles between various logarithmic terms. Here, 
we do not want to disentangle the implicit consistency conditions. Rather, with an 
understanding of which spurious poles appear in which integral coefficient, we try to 
maneuver around them. A detailed discussion of the factorisation properties of one-loop 
amplitudes, as well the spurious singularities that appear, may be found in refs. [22, 140]. 

The factorisation of amplitudes follows from the combined behaviour of integral 
functions and the integral coefficients in the factorisation limit. If wc turn this around, 
given the general factorisation (102) of an amplitude and given factorization properties 
of integral functions [22, 140], we may then determine the factorisation properties of 
the integral coefficients. By applying this logic to multi-particle factorisations (102) 
we conclude that the coefficients behave as if they were tree amplitudes as long as 
the factorisations are entirely within a cluster of legs (and are not on the momentum 
invariant of the entire cluster). That is the coefficient behaves as. 



with notation as indicated in fig. 12. For convenience, from now on we will use the 
notation q for all integral coefficients; coefficients of two-point functions, three-point 
functions etc. 

Assuming that the spurious denominators do not pick up a z dependence — below 
we describe simple criteria for ensuring this — we obtain a recursion relation for the 




(109) 



h 
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Figure 12. Diagrammatic contributions to on-shell recursion for a three-point integral 
coefBcient Ci_„{z) for a [j, I) shift. The left-hand-side indicates a triangle loop diagram 
that gives rise to the input lower-point coefficient Ci^rn+i- The recursion is built from 
sewing a lower-point integral coefficient Ci^rn+i{za) and a tree amplitude An-m+iiza) 
The vertices of the triangle graph are shown as the tree amplitudes of the underlying 
triple cut. The shifted legs are picked from within the same cluster of legs. 



coefficients which strikingly is no more comphcated than for tree amphtudes, 

a,h " 

where ^n-m„ -i-i(^a) Ci m^+iiza) are shifted tree amphtudes and coefficients evaluated 
at the residue value , h denotes the hehcity of the intermediate state corresponding 
to the propagator term i / . In this expression one should only sum over a limited 
set of poles; if the shifts are chosen from within a cluster, the only poles that should be 
included are from within the kinematic invariants formed from the momenta making up 
the cluster. Pictorially, this coefficient recursion relation is shown in fig. 12. 

Several consistency requirements have to be fulffiled for on-shell recursions of 
integral coefficients. Simple criteria for a valid recursion were given in [20]: 

(i) The shifted tree amplitude, on the side of the cluster undergoing recursion, vanishes 
for large z. 

(ii) The loop- momentum dependent kinematic poles are unmodified by the shift (96). 

In addition to the standard requirement, that the shifted coefficient vanishes for large 
z, we have criteria (ii) which sidesteps contributions from spurious poles. In particular, 
this relies on the assumption that spurious poles in z appear whenever loop-momentum 
dependent propagator give rise to non-zero residues. 

The above criteria are fulfilled for selecting a shift within the same cluster of legs 
in addition to a requirement on the helicity structure of the respective cluster. For 
a particular set of helicity amplitudes, so called split helicity amplitudes, coefficient 
recursions take a very simple form. (Split helicity refers to color-ordered amplitudes with 
all like helicities adjacent; y4„(— . . . — , -|- . . . -|-).) For these amplitudes computations 
were performed for all multiplicity computations in supersymmetric Yang-Mills, pure 
QCD as well as all 7-point amplitudes for W/Z + 3-}et production [30, 32, 34], including 
one as well as two fermion lines. Further examples include integral coefficients with 
split-helicity corners. 
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parent : daughters 




Figure 13. Parent and daughter integral functions. The Gram determinant 
singularity associated to the parent 3-point integral is A(Ki2, k^, K^c^). It is present in 
its integral coefficient as well as in the coefficients of both daughter 2-point integrals. 



5.5. Spurious Recursions. 

When recursions for integral coefficients can be setup, we can formulate analogous 
recursions for rational terms. These 'auxiliary recursions' then directly give access 
to spurious residues within rational terms. We obtain in this way a purely recursive 
approach for rational terms, where no reference to integral expansions of cut pieces is 
needed. 

In fact, given that we shift integral coefficients within one of its corners, integral 
functions stay inert under the shift; they only depend on un-shifted momenta. Similarly, 
we may expand integrals around their Gram-determinant poles (see eq. (107)), still 
yielding terms independent of the shift-variable. Similar steps can be followed for 
integral coefficients of lower-point (daughter) integrals which inherit a given Gram 
determinant singularity. The structure of parent and daughter integral coefficients is 
indicated in fig. 13. 

Putting all pieces together, we obtain an on-shell recursion relation for rational 
terms, which is valid close to a given vanishing Gram determinant, 

i i,a,h 

where Ajs denotes the Gram determinant around which we intend to expand the rational 
term. Xf denote the integral functions, which are expanded around the location of the 
zero of the Gram determinant as in eq. (107). In addition, only their rational terms are 
kept. The symbol a labels the physical poles that need to be considered in the coefficient 
recursions of the parent integral coefficient. For convenience integral coefficients are 
denoted by q, for any of the integral functions somewhat differing from our earlier 
notation in section 4. 

We thus obtain a recursion relation for the rational terms, which is exact near the 
Gram determinant singularities. Contour integrals whether continuous or discrete can 
then be used to extract the exact value of the spurious residues of eq. (105). 

Thus obtained recursive approaches allow to compute rational terms from purely 
recursive methods, albeit, using a set of well chosen auxiliary shifts. These methods have 
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already been applied for the computations of split helicity scattering amplitudes, needed 
in the computation of W^ + 3-jet and Z + 3-jet production [30, 32, 34] leading to greatly 
improved computation time for these pieces. A much wider class of computations seems 
in reach using the above newly developed recursive methods. 



6. Conclusions 

The computations of observable cross-sections rely on the composition of several 
components, which are linked through the fundamental QCD Lagrangian. Here we 
described modern approaches for evaluating the hard scattering part at ncxt-to-lcading 
order, which can be programmed. The value of a first principle understanding of 
scattering processes in addition to the increased quantitative control motivates the quest 
for cross sections at NLO. 

The description of processes with complex final states at NLO is one of the central 
achievements of many recent developments of numerical on-shcll [47, 19, 20, 21] and 
unitarity methods [22, 23, 25, 27, 21, 29]. These methods are already used by a 
new generation of tools [79, 21, 68, 133, 39, 40, 41, 38, 42]. Such numerical methods 
scale very well as the number of external states increases and have already lead to a 
new understanding for several proton scattering processes at hadron colliders [33, 35, 
36, 37, 38] including W/Z + 3-jet and W + 4-jet production [30, 31, 32, 36]. Being 
key backgrounds to many new physics signals including supersymmetry searches, the 
explicit results emphasizes the importance of modern field-theory methods for precision 
prediction of hadron collider physics. 

Numerical on-shell and unitarity methods allow us to automate computations of 
scattering amplitudes in a numerically stable and efficient way. These methods exploit, 
for example, the discontinuities across branch cuts to construct amplitudes. Efficiency 
and stability originates then in two facts: Firstly, discontinuities are expressed in terms 
of purely on-shell information, namely, on-shell tree- amplitudes. Evaluation of on-shell 
trees is fast and numerically stable. In addition, this allows us to ignore ghosts and 
gauge fixing altogether and reduces the use of redundant unphysical information. 

These modern ideas presented here my also help with issues beyond the computation 
of matrix elements. New insights may impact on various other components of multi- 
jet computations. This includes subtraction methods for integrating parton level NLO 
computations, showering and formulation of new observables. 

The central strategy, leading to the success and popularity of these developments, 
is to make maximal use of physical principles and mathematical structures in order to 
obtain efficient and robust phenomenology. This kind of work, is at the cross-roads 
of theory and phenomenology has led to many recent insights, on the more formal 
side [146, 147, 21, 29, 136, 148] as well as of more phenomenological nature [110]. 
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